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ABSTRACT 


A  review  of  basic  pulsar  observations,  complete  until  early 
1979,  is  presented,  followed  by  the  arguments  leading  to  the 
identification  of  neutron  stars  as  the  seat  of  the  pulsar  phenomenon. 

The  basic  properties  of  neutron  stars  as  deduced  from  the  Tolman- 
Oppenheimer-Volkof f  equation  and  current  ideas  of  high-density 
equations  of  state  follow.  The  structure  of  matter  at  the  surface  of 
a  magnetized  neutron  star  is  presented  in  a  clear  and  consistent 
manner,  and  a  mechanism  to  prevent  the  ionization  of  free  atoms  on 
the  surface  is  proposed.  Due  to  high  binding  energies  and  this 
mechanism  it  is  shown  that  the  neutron  star  surface  must  be  ion- 
retentive  over  most  of  the  star.  The  equations  governing  inertia-free 
plasma  around  the  neutron  star  are  presented  and  the  vacuum  electric 
field  is  found.  The  plasma  charge  density  for  a  force-free  magneto¬ 
sphere  and  its  physical  conditions  are  discussed.  Possible  departures 
from  force-free  conditions  and  effects  of  stellar  charge  and  surface  ion 
retention  are  noted,  and  the  formation  of  polar  gaps  above  an  aligned 
rotator  where  wy  <  0  and  charge  is  swept  from  polar  regions  is  shown 
to  be  caused  by  the  ion  retention,  according  to  the  Ruderman-Sutherland 
theory.  The  characteristics  of  discharges  due  to  pair  production  are 
discussed.  It  is  proposed  that  synchrotron  radiation  from  the  produced 
pairs,  and  annihilation  radiation  from  ejected  positrons,  will  be 
observed  from  pulsars  when  observational  techniques  permit.  Radiation 
mechanisms,  especially  those  due  to  pair  production,  are  briefly 
discussed. 
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CHAPTER  1 


BASIC  PROPERTIES  OF  PULSARS  -  OBSERVATIONS  AND  MODELS 

Pulsars  have  aroused  a  large  amount  of  interest  since  their 
discovery  twelve  years  ago  (Hewish  et  al.  1968).  As  a  result  there 
exists  an  extensive  set  of  detailed  observations  of  them,  backed  by  a 
bewildering  array  of  papers  on  pulsar  theories.  Despite  this  huge 
amount  of  effort  in  both  observation  and  explanation,  no  fully 
satisfactory  theory  or  model  has  yet  been  made  to  explain  the  details 
of  pulsar  emission.  Indeed  there  has  been  a  tendency  to  construct 
models  independently  explaining  single  aspects  of  the  phenomenon,  with 
hopes  that  these  narrower  views  are  not  contradictory. 

In  making  pulsar  models,  one  uses  established  physical  theory, 
guided  by  observational  constraints,  to  arrive  at  a  consistent  solution. 
This  may  in  turn  predict  further  properties,  which,  if  observed,  lend 
more  weight  to  the  model.  Since  present-day  models  are  limited  in 
detail,  we  should  examine  the  general  properties  which  must  be 
explained  if  we  wish  to  claim  any  understanding  at  all  of  the  pulsar 
mechanism. 

More  detailed  discussions  of  observations  are  presented  in  two 
recent  books  (Manchester  and  Taylor  1977;  Smith  1977). 

Pulses 

Regular  pulses  of  enhanced  emission  are  the  defining  character¬ 
istic  of  pulsars.  The  emission  is  from  objects  which  are  starlike  in 
angular  size,  both  at  optical  frequencies  (see  e.g.  Peterson  et  al. 
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1978),  and  in  the  radio  region  after  allowance  for  interstellar 
scintillation  (Smith  1977;  Andrew  et  al.  1964;  Coles  and  Kaufman  1977). 
These  pulses  have  broadband  spectra  and  for  various  pulsars  are 
observed  in  spectral  regions  from  radio  to  gamma-ray.  The  pulses  are 
themselves  composed  of  one  or  more  subpulses  which  occur  with 
irregular  timing  and  intensity,  but  such  that  a  time-averaged  intensity 
profile  of  a  large  number  of  pulses  maintains  an  extremely  stable  shape, 
although  one  unique  to  each  pulsar.  With  higher  time  resolution,  the 
subpulses  of  5  of  8  pulsars  so  studied  (Cordes  1975)  showed  micro¬ 
pulses  on  a  time  scale  of  hundreds  of  microseconds.  Secular  changes  in 
this  stable  integrated  profile  have  been  seen  only  in  the  case  of  the 
single  known  binary  radio  pulsar  (Taylor  et  al.  1978).  The  period  of 
the  integrated  profile  is  extremely  stable,  to  as  much  as  one  part  in 
10 12  over  several  years.  Periods  ranging  from  .033  s  to  4.308  s  have 
been  observed.  First  derivatives  of  the  period  range  from  +10~18ss-1 
to  +10~12ss_1  for  radio  pulsars,  and  are  larger  in  magnitude  but 
negative  in  sign  for  the  X-ray  (binary)  pulsars.  Abrupt  period 
decreases  known  as  "glitches”  have  been  observed  in  those  two  pulsars 
considered  to  be  the  youngest,  the  Crab  pulsar  P0531+21,  and  the  Vela 
pulsar  P0833-45,  and  also  in  P1641-45  (Manchester  1979).  Small 
glitches’  in  a  loose  sense"  have  been  observed  in  nine  other  pulsars 
in  very  accurate  pulse  timing  studies  (Gullahorn  and  Rankin  1978b) . 

The  glitches  of  the  Vela  pulsar,  illustrated  in  Figure  1,  have 
been  four  in  number  (the  last  in  July  1978),  each  a  period  decrease  of 
about  200  nanoseconds,  a  fractional  change  of  2  x  10“6  in  period. 
Glitches  of  the  other  two  pulsars  have  been  a  factor  of  100  smaller, 
and  in  the  Crab  vary  in  size.  The  period  derivative  decreases  during 
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Figure  1.  Glitches  of  the  Vela  Pulsar  super¬ 
posed  on  the  continuously  increasing  period, 
(adapted  from  Manchester  and  Taylor  1977 , 
p.  114,  and  Manchester  1979). 


a  glitch,  the  fractional  change  being  measured  in  tenths  of  a  percent. 
The  period  derivative  returns  to  its  pre-glitch  value  on  a  time  scale 
of  several  days  in  the  Crab  pulsar,  and  of  about  a  year  in  the  Vela 
pulsar.  Higher  derivatives  and  irregular  timing  effects  have  been 
noted  in  most  of  the  better-observed  pulsars. 

In  most  pulsars  the  radio  emission  occurs  during  only  about 
3  percent  of  the  cycle,  but  some  pulsars  exhibit  an  "interpulse" 
roughly  in  antiphase  with  the  main  pulse  (Figure  2) .  "Marching 
subpulses,"  the  correlation  of  the  position  of  a  subpulse  in  one  cycle 
to  its  position  in  the  next  cycle  (Figure  3),  and  "mode  changing,"  in 
which  two  different,  stable  integrated  profiles  are  observed  at 
different  times,  characterize  a  small  number  of  pulsars.  At  optical 
frequencies,  unpulsed  emission  has  been  observed  from  the  Vela  and 
Crab  pulsars  (Peterson  et  al.  1978).  At  radio  frequencies  the  unpulsed 
component  is  small  (Huguenin  et  al.  1971). 

Although  the  integrated  profiles  are  stable  in  shape,  their 
intensity  may  vary  considerably  on  time  scales  of  hours  or  more,  in 
some  cases  going  below  detectable  levels  for  long  periods  of  time. 

The  pulsed  emission  usually  has  very  high  linear  polarization 
whose  position  angle  changes  mono tonically  across  the  integrated  pulse. 
This  polarization  decreases  with  increasing  frequency. 

Pulsar  emission  peaks  in  the  radio  region  and  there  the  flux 
density  S  has  a  power  law  spectrum  of  form: 

S(v)  =  SQva  (1-1) 

where  SG  =  S(v0)/v0a,  v  is  frequency,  vG  a  reference  frequency,  and 
the  index  a  typically  about  -1.5.  (For  a  thermal  source,  the  index  a 


5 


Figure  2.  Pulse  profiles  of  the  Crab  Pulsar 
P0531+21  at  optical  (above)  and  430  MHz  (below) 
frequencies.  Very  short  timescale  structures 
at  radio  frequency  may  be  seen.  (Drake  1970) . 
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Figure  3.  Drifting  subpulses.  The  lines  of 
horizontal  dots  represent  individual  pulses 
the  dot  size  reflects  the  intensity  at  a  given 
longitude  in  the  pulse.  The  integrated  profiles 
are  shown  at  the  top.  (Manchester  1979). 


would  have  value  +2 . ) 


The  binary  X-ray  pulsars,  whose  spectra  peak  sharply  in  the 
X-ray  region,  do  not  follow  this  power  law  spectrum,  but  the  fact  that 
at  high  frequencies,  only  the  Crab,  Vela  and  two  longer-period  pulsars 
have  been  observed,  is  consistent  with  this  steep  spectrum  for  all 
radio  pulsars.  (Optical,  X-ray,  and  y-ray  pulses  have  been  detected 
from  the  Crab  pulsar,  P0531+21,  optical  and  y-ray  pulses  from  the  Vela 
pulsar  P0833-45,  and  y-ray  pulses  from  P1747-46  and  P1818-04 
(Ogleman  et  al.  1976;  Buccheri  et  al .  1979).) 

In  the  binary  X-ray  pulsar  Her  X-l,  a  strong  emission  line  at 
53  keV  has  been  observed  superposed  upon  an  X-ray  exponential  spectrum 
(Trumper  et  al.  1977;  Levi  1977).  The  intensity  of  the  line  is  such 
that  atomic  or  nuclear  sources  seem  very  unlikely,  and  it  has  been 
interpreted  as  electron  cyclotron  ("non-relativistic  synchrotron") 
radiation  from  electrons  moving  in  a  magnetic  field  of  strength 
5  x  1012  G. 


The  high-frequency  synchrotron  radiation  emitted  by  the  Crab 
Nebula  is  steady  in  intensity  for  many  years  despite  the  fact  that  the 
electrons  producing  it  have  much  shorter  emitting  lifetimes.  The 
continuous  injection  of  1038  erg*s-1  of  energy  carried  by  relativistic 
electrons  from  the  Crab  pulsar  is  the  only  plausible  solution.  The 
expansion  of  the  nebula  accelerates  at  a  rate  implying  that  a  further 
4  x  1038  erg*s_1  are  pumped  into  gas  kinetic  energy.  The  total 
luminosity  of  the  pulsar  must  thus  be  at  least  5  x  1038  erg*s_1. 


;  f  IK  Si 


The  luminosity  of  pulsar  emission  cannot  be  directly  measured 
because  assuming  isotropic  emission  (as  is  usually  done  with  stars) 
does  not  seem  justifiable.  Thus,  although  flux  can  be  measured,  and  the 
fact  that  the  speed  of  radio  waves  in  the  interstellar  medium  is 
frequency-dependent  allows  distances  to  be  assigned,  the  luminosity  of 
a  pulsar  depends  upon  the  model  chosen  and  only  in  this  sense  is  it 
"observable."  Thus  only  the  Crab  Nebula  and  the  possibility  that  high 
energy  cosmic  ray  electrons  are  pulsar-produced  (Ostriker  1972)  allow 
luminosities  to  be  determined.  More  data  on  pulsar  luminosity,  from 
sources  other  than  the  Crab  Nebula,  would  be  very  desirable. 

The  Crab  Nebula  contains  a  widespread  magnetic  field  of  one 
milligauss  according  to  the  interpretation  of  its  continuum  radiation 
as  synchrotron  radiation  from  electrons.  This  nebular  field  is  too 
large  in  magnitude  and  extent  to  be  explained  without  invoking  the  pulsar 
as  its  source. 

Of  limited  physical  significance  without  a  model,  but  useful 
statistically  in  estimating  the  degree  of  completeness  in  pulsar 
surveys,  is  the  luminosity  function  of  pulsars  in  the  solar  neighbor¬ 
hood.  Defining  the  luminosity,  L,  as 

L  =  S4oo<i2  (1~2) 

where  S4qo  is  the  400  MHz  flux  density  (in  Jy)  and  d  the  distance  in 
parsecs,  the  z-integrated  luminosity  function  per  semidecade  luminosity 
interval. is  (Manchester  and  Taylor  1977,  p.  150): 

*  1 2  kpc~2 


$(L)  =  350L 


(1-3) 


Mechanical  Properties 


A  study  of  the  binary  pulsar  P1913+16  (Taylor  et  al.  1978)  has 
revealed  effects  predicted  by  General  Relativity  (GR)  which  could, 
however,  instead  be  explained  as  tidal  effects  in  the  pulsar  or  its 
companion.  If  these  effects  are  indeed  from  GR,  a  very  rigid  structure 
must  characterize  the  two  objects. 

A  mass  of  1.39  ±  .15  M^  for  the  pulsar  and  1.44  ±  .15  M0  for  the 
companion  were  also  determined. 

In  an  analysis  of  5  binary  X-ray  pulsars  whose  companions  could 
be  observed  optically  (Rappaport  and  Joss  1977b) ,  pulsar  masses  ranging 
from  0  to  4  Mg  in  individual  cases  are  not  observationally  excluded. 
However,  a  common  mass  of  1.2  to  1.8  M@  would  be  consistent  with 
observations  of  all  of  the  systems. 

Statistical  Properties 

As  of  late  1978,  321  pulsars  had  been  observed,  241  of  them  in 
surveys  specifically  designed  to  provide  a  meaningfully  uniform 
statistical  sample  (Manchester  et  al.  1978;  Damashek  and  Taylor  1978). 
Sources  for  statistical  data  presently  available  are  appended  to  this 
thesis. 

The  observed  distribution  of  periods  is  shown  in  Figure  4,  and 
since  searches  have  covered  fairly  equally  periods  in  the  range  .01  to 
10  seconds,  which  is  plotted,  this  is  likely  to  represent  the  true 
distribution  of  pulsar  periods.  The  distribution  is  apparently 
unimodal:  the  slight  dip  at  P  =  1  s,  which  in  less  complete  data 

appeared  to  indicate  bimodality  (Manchester  and  Taylor  1977,  pp.  152-153) 
no  longer  seems  significant. 


NUMBER  OF  PULSARS 


Figure  4.  Distribution  of  periods  among 
321  radio  pulsars.  (Manchester  1979) 


Period  derivatives,  when  measured,  are  always  positive  for 
radio  pulsars  -  their  period  steadily  and  slowly  lengthens.  As  Figure  5 
indicates,  there  is  no  pronounced  correlation  between  period  and 
derivative.  (Note  that  the  graph  is  logarithmic  on  both  axes.) 

An  important  selection  effect  which  enters  this  graph  is  that 
only  pulsars  with  very  large  derivatives  or  those  which  have  been 
carefully  observed  for  sufficiently  large  amounts  of  time  have  had  their 
derivatives  observed.  The  smaller  derivatives  are  thus  discriminated 
against. 

A  characteristic  age  x  =  ^P/P  may  be  defined  and  be  expected  to 
be  an  order  of  magnitude  estimate  of  the  true  age  of  a  pulsar.  The 
distribution  of  ages  shown  in  Figure  6  may  have  larger  ages 
discriminated  against  due  to  the  difficulty  of  observing  small 
derivatives. 

Figure  7a  clearly  shows  the  clustering  of  pulsars  toward  the 
galactic  equator,  proving  that  they  are  generally  in  the  disk  of  our 
galaxy.  Figure  7b,  showing  pulsars  more  than  one  kiloparsec  from  the 
Sun  projected  on  the  galactic  plane,  shows  that  surveys  beyond  about 
5  kpc  cannot  be  very  complete,  but  that  there  is  a  real  tendency  for 
the  pulsar  density  distribution  to  peak  inside  the  galactic  radius  of 
the  Sun.  Figure  7,  from  Manchester  (1979),  is  based  on  pulsar 
distances,  d,  calculated  from  dispersion  measures  (DM) 

DM  =  /n  dl  (1-4) 

o  e 

which  can  be  measured  due  to  the  delay  in  pulse  arrival  times  at 
frequencies  mi  and  002 »  due  to  interstellar  dispersion: 


Period  Derivative  (s  s 
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Figure  5.  Period  derivative  versus  period  for 
103  pulsars.  Solid  lines  have  constant  PP, 
dashed  lines  constant  characteristic  age. 

Solid  dots  are  data  from  Manchester  and  Taylor 
1977,  p.  110,  empty  dots  data  from  Gullahom 
and  Rankin  1978. 
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Figure  6.  Distribution  of  83  pulsar  character¬ 
istic  ages,  T  ,  less  than  2  x  10 ^  a.  (24  others 
have  T  >  2  x  10^  a) .  Data  are  from  Manchester 
1979,  Gullahom  and  Rankin  1978b. 
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90° 


Figure  7.  (a)  Distribution 
in  1  and  b  of  321  pulsars. 


180° 


(b)  Distribution  in  1 
and  R  of  pulsars  with 
d>  1  kpc.  (Manchester 
1979) 
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?  2 

t2  “  ti  =  (a>2_2  -  o)1_2)DM  (1-5) 

With  a  model  for  the  electron  density  distribution,  ne  (r,l,b),  such 
as  the  one  used  in  Manchester  (1979),  consisting  of  a  uniform  0.025  cm-3 
superposed  on  an  exponential  layer  0.015e  cm-3  and  known  local 

HII  regions,  the  distance  may  be  uniquely  determined  from  the  DM. 

The  distribution  in  z  shown  in  Figure  8  can  be  fitted  by  an 
exponential  with  scale  height  400  pc.  This  is  much  larger  than  that 
of  0-B  stars  (80  pc)  or  SNR  (60  pc).  Pulsars  are,  however,  high- 
velocity  objects  with  an  "unbiased”  mean  z-velocity  of  60  km  s-1 
(Manchester  1979)  based  on  measured  proper  motions  (Manchester  and 
Taylor  1977;  Gullahorn  and  Rankin  1978a).  At  this  velocity,  the  mean 
z  of  300  pc  would  be  reached  (from  z  =  0)  in  about  5.5  x  106  a,  which 
is  of,  the  order  of  the  mean  characteristic  age  of  the  pulsars  shown  in 
Figure  8. 

Among  radio  pulsars  only  P1913+16  among  321  detected  is  a  known 
member  of  a  binary  system,  although  for  about  half  of  these  detections 
the  data  are  not  yet  complete  enough  to  exclude  some  being  binaries 
(Manchester  et  al.  1978).  At  least  half  of  the  "normal"  stars  of  our 
Galaxy  are  in  binary  systems. 

Basic  Characteristics  of  Pulsar  Models 

The  most  striking  feature  of  the  observational  data  is  the  very 
fact  that  pulsed  emission  from  astronomical  objects  can  occur  on  a  very 
short  time  scale  and  with  great  regularity.  The  short  period  is  so 
unique  in  astronomical  contexts  that  an  argument  can  be  made  that  the 
reason  that  pulsars  were  not  discovered  until  1967  is  that  no  one 
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Figure  8.  z-distribution  of  321  pulsars. 
(Manchester  1979) 


looked.  Indeed  P0329+54  and  the  Crab  pulsar  were  both  recorded  "before" 
pulsars  were  discovered  but  their  unique  signals  not  recognized  (Smith 
1977).  After  Jocelyn  Bell’s  discovery  in  July  1967  of  the  signals  of 
P1919+21  on  recordings  taken  for  the  study  of  interstellar 
scintillation,  the  group  at  Cambridge,  of  which  she  was  part,  at  first 
hypothesized  that  their  "remarkable  nature"  was  due  to  their  generation 
by  man  (Hewish  et  al.  1968)  or  by  "some  silly  lot  of  little  green  men" 
(Burnell  1977).  However,  the  absence  of  parallax,  code,  or  source’s 
orbital  Doppler  shift,  and  the  subsequent  discovery  of  three  further 
pulsars  (by  the  Cambridge  group  in  the  cautious  seven  months  between 
the  initial  discovery  and  publication) ,  indicated  that  the  phenomenon 
was  natural.  The  pulse  widths,  of  order  10  ms,  led  on  the  basis  of 
light  travel-time  arguments  to  the  conclusion  that  the  sources  would  be 
of  planetary  or  smaller  size,  and  the  discovery  announcement  (Hewish  et 
al.  1968)  itself  proposed  that  white  dwarfs  or  neutron  stars  generated 
the  emission. 

Pulsation  of  white  dwarfs  was  initially  favoured  as  the  basic 
timing  mechanism,  pulsation  periods  as  short  as  .07  s  being  possible 
for  non-uniformly  rotating  white  dwarfs  (Durney  et  al.  1968).  The 
discovery  of  the  Crab  and  Vela  pulsars,  with  considerably  shorter 
periods,  ruled  out  white  dwarfs,  and  as  neutron  star  pulsation  periods 
are  far  too  short,  pulsation  became  disfavoured  as  a  timing  mechanism 
(Manchester  and  Taylor  1977).  Theories  invoking  pulsation  or  explosive 
radial  motion  of  only  surface  layers  of  white  dwarfs  or  neutron  stars 
(Hoyle  and  Narlikar  1968;  Israel  1968a,  b)  also  predicted  periods  of 


about  one  second. 


Orbital  motion  about  a  neutron  star  could  occur  with  periods 


down  to  about  1  ms  and  possibly  cause  the  pulses  by  gravitational 
focusing  (Saslaw  et  al.  1968)  but  gravitational  radiation  would  cause 
the  period  to  decrease,  and  the  tidal  forces  acting  on  a  satellite 
would  require  it  to  have  unusual  properties  to  remain  intact  (Smith 
1977)  .  The  discovery  of  increasing  periods  for  the  Crab  pulsar 
(Richards  and  Cornelia  1969),  and  subsequently  for  all  others,  led  the 
orbital  theory  to  disfavour. 

The  remaining  mechanism  involving  condensed  objects  is  that  the 
rotation  of  a  massive  body  acts  as  the  timing  standard,  just  as  (until 
recently)  on  Earth.  White  dwarfs  are  stable  against  centrifugal 
breakup  only  for  periods  longer  than  one  second  (Ostriker  1968)  and  so 
were  ruled  out.  Pacini  (1967)  had  suggested  that  a  rotating, 
magnetized  neutron  star  could  emit  large  amounts  of  energy  in  the  form 
of  low-frequency  electromagnetic  waves,  and  Gold  (1968)  was  the  first 
to  propose  such  a  star  as  the  timing  source  in  pulsars,  and  correctly 
predict  that  positive  period  derivatives  should  be  seen  due  to  the 
slowing  of  the  rotation  as  energy  radiated  from  the  star. 

By  the  middle  of  1969  there  was  general  agreement  that  Gold’s 
model  provided  a  basic  timing  mechanism  with  all  of  the  needed 
properties.  In  the  intervening  ten  years  the  problems  of  understanding 
the  atmosphere  of  the  neutron  star  and  pulse  production  have  "proven 
to  be  much  more  difficult  than  might  have  reasonably  been  anticipated" 
(Jackson  1976).  The  remaining  chapters  of  this  thesis  address  these 
problems. 

For  completeness,  some  post-1969  models,  which  are  considered 
"non-standard"  as  opposed  to  the  rotating  neutron  star  hypothesis,  will 


be  briefly  discussed. 


Similar  in  some  respects  to  pulsational  theories  is  the  "bouncing- 
core  theory"  (Harrison  1970),  in  which  the  core  of  an  object  on  the  verge 
of  gravitational  collapse  is  unstable  and  does  collapse  periodically 
until  a  harder  equation  of  state  in  the  collapsed  state,  and  changes 
in  the  metric  around  the  object,  halt  and  reverse  the  collapse.  This 
model  might  allow  periodicities  of  the  right  order,  but  the  energy 
radiated  through  the  non-collapsing  crust  (which  is  supported  by 
radiation  pressure)  would  likely  be  more  like  that  of  a  blackbody  than 
that  observed  in  pulsars. 

The  solar  cycle  theory  advanced  by  Apparao  and  Chitre  (1970) 
attributes  pulsar  emission  to  activity  similar  to  solar  flares, 
associated  with  a  neutron  star  magnetic  cycle  similar  to  that  of  the 
Sun  but  on  a  much  shorter  time  scale.  The  precision  of  the  periodicity 
and  time  derivatives  are  not  convincingly  explained  by  such  a  theory. 

The  mode-locked  maser  theory  of  Mertz  (1974)  explains  pulse 
shapes  and  periods  through  a  coherent  mechanism  in  the  resonant  cavity 
formed  by  the  surface  of  a  white  dwarf  star  and  its  hypothesized  shell¬ 
like  ionosphere.  The  supply  of  energy  to  such  a  resonant  system  for 
coherent  emission  is  this  theory’s  major  difficulty. 

The  major  gaps  in  these  theories  lend  support  to  the  neutron 
star  hypothesis  through  "what-else-can-it-be"  arguments.  The  ability 
of  a  rotating  neutron  star  to  explain  timing  observations,  supply 
energy  at  appropriate  rates,  and  explain  certain  radiation  features, 
has  led  to  the  conviction  that  further  work  un  it  will  eventually  bear 
fruit,  even  though  present  models  based  on  it  are  by  no  means  fully 
sat:  ^factory. 


CHAPTER  2 


STRUCTURE  OF  NEUTRON  STARS 


The  aim  of  this  brief  discussion  is  to  show  that  neutron  stars 
could  be  expected  to  have  masses  and  moments  of  inertia  consistent  with 
their  being  pulsars,  and  to  discuss  magnetic  fields  and  boundary 
conditions  which  will  be  needed  for  the  formulation  of  the  pulsar 
mechanism  as  an  electrodynamical  problem. 


Neutron  Stars 

Very  soon  after  the  discovery  of  the  massive  but  electrically 
neutral  neutron  it  was  realised  that  a  large  number  of  neutrons  could  be 
held  together  gravitationally,  with  pressure  support  essentially  due  to 
the  repulsive  core  of  the  nucleon  potential,  to  form  a  stable  condensed 
object.  Oppenheimer  and  Volkoff  in  1939  made  the  first  neutron  star 
model,  and  since  the  discovery  of  pulsars  in  1967  and  their 
"identification"  as  rotating  neutron  stars  in  1968  (Gold  1968)  much 
work  has  been  done,  to  the  point  that  most  models  agree  on  the  general 
structure  and  differ  only  in  detail,  and  thus  may  be  considered 
satisfactory  (Canuto  1977). 

Within  the  framework  of  General  Relativity  (GR)  the  equation  of 
hydrostatic  equilibrium  for  a  non-rotating  spherically  symmetric  neutron 
star  is  the  Tolman-Oppenheimer-Volkof f  (TOY)  equation 


dP 

dr 
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where  P(r)  is  the  pressure,  and  p(r)  the  density,  at  radius  r, 

m(r)  =  /^4ttp  (r  ’  )r '  2dr  ’  (2-2) 

being  the  mass  interior  to  r. 

Other  gravitational  theories  than  GR  would  cause  differences 
from  TOV  solutions  which  vary  from  being  very  minor  (Brans-Dicke  theory) 
to  non-negligible  or  even  very  large  (Brecher  and  Caporaso  1977) . 
However,  such  is  the  present  status  of  GR  that  the  use  of  the  TOV 
equation  has  formed  the  core  of  recent  neutron  star  theory  (Canuto 
1977) .  The  consistency  of  neutron  star  and  pulsar  parameters  with  TOV 
predictions  and  with  GR  will  be  examined  later  in  this  chapter. 

As  is  the  case  in  non-relativistic  stellar  modelling,  the  use  of 
an  equation  of  state 


P  =  P(p)  (2-3) 

and  a  boundary  condition  P(R)  =  0  determines  a  unique  total  mass 
M  =  m(R) ,  and  total  radius  R  for  a  chosen  central  density  pc. 

Five  regions  of  distinctly  different  properties  and  equations 
of  state  may  be  identified  in  terms  of  density  regimes  (Lamb  1977) : 

1.  The  surface  (p  <  106  g  cm”3)  where  the  effects  of 

temperature  and  electromagnetic  fields  can  be 
significant. 

2.  The  outer  crust  (106  g  cm”3  <  p  <  4  x  1011  g  cm”3); 

Here  a  solid  lattice  of  nuclei  is  in  a  sea  of 
degenerate  relativistic  electrons,  as  in  white  dwarf 
interiors.  Iron  (56Fe),  the  most  stable  nucleus  at  low 


pressure,  will  be  the  most  abundant,  until  above 
8  x  10 6  g  cm-3  more  neutron  abundant  nuclei 
predominate  due  to  their  inability  to  decay  by 
releasing  electrons  into  a  full  Fermi  sea.  At  about 
4  x  1011  g  cm-3  is  the  "neutron  drip  point"  where 
the  nuclei  are  neutron  saturated  and  free  neutrons 
exist. 

3.  The  inner  crust  (4  x  1011  g  cm-3  <  p  <  2  x  1014  g  cm"3), 
below  the  drip  point,  contains  a  proportion  of  neutrons 
increasing  rapidly  with  pressure  until  the  lattice 
dissolves  and  a  sea  of  neutrons,  with  only  a  few  percent 
free  electrons  and  protons,  exists. 

4.  The  neutron  liquid  exists  above  p  =  2  x  1014  g  cm-3 
until  some  possible  critical  density. 

5.  A  core  may  exist,  consisting  of  solid  neutrons,  hyperon 
or  quark  "soup,"  or  pion  condensate  (Baym  and  Pe thick 
1975),  for  p  >  1015  g  cm-3  (Manchester  and  Taylor  1977). 

Integration  of  the  TOV  equation  uses  equations  of  state, 

P  =  P(p),  appropriate  to  these  various  regimes.  The  pressure  in  a 
neutron  star  is  almost  entirely  due  to  zero-point  energy,  and  thus  any 
possible  temperature  dependence  is  ignored  (Smith  1977,  p.  36).  Possible 
effects  of  temperature  at  the  surface  will  be  discussed  later. 

A  recent  review  (Canuto  1977)  of  10  TOV-based  neutron  star 
models  showed  that  there  is  agreement  among  them  in  predicting  masses 
of  up  to  2  M@,  moments  of  inertia  in  the  range  1044-1045  g  cm2,  and 
radii  of  8  to  12  km.  These  masses  would  be  consistent  with  those  of 


pulsars  discussed  in  Chapter  1. 

From  the  understanding  of  the  properties  of  matter  in  these 
five  regions  come  predictions  of  further  mechanical  properties  of 
neutron  stars.  The  iron  lattices  at  the  surface  and  in  the  outer  crust 
will  be  very  rigid  due  to  the  large  increases  in  Coulomb  energy 
associated  with  any  deformation  from  the  lowest  energy  position. 

The  neutrons  in  the  inner  crust  and  deeper  are  superfluid  (Manchester 
and  Taylor  1977),  and  are  very  weakly  coupled  to  the  crust. 

The  result  of  the  rigid  outer  structure  and  superfluid  interior 
is  that  "crustquakes"  occur.  As  the  neutron  star  loses  rotational 
energy,  its  equilibrium  shape  changes.  However,  the  rigidity  of  the 
outer  crust  does  not  allow  it  to  deform  continuously  -  instead  stresses 
build  up  until  a  sudden  deformation  occurs,  resulting  in  a  decrease  in 
(rotation- induced)  oblateness  and  in  moment  of  inertia.  The 
conservation  of  angular  momentum  then  forces  a  decrease  in  period  of 
the  crust.  Due  to  the  small  coupling  torque  between  crust  and  super¬ 
fluid,  which  is  proportional  to  their  relative  angular  velocity,  the 
torque  slowing  the  crust  is  increased,  and  thus  so  is  the  rotation 
frequency  derivative  -  the  period  derivative  decreases,  as  is  observed. 
As  the  interior  and  crust  come  to  equilibrium,  the  torque  vanishes 
and  the  period  derivative  returns  to  its  previous  value.  This  theory 
predicts  the  order  of  magnitude  of  the  time  interval  between  glitches 
in  the  Crab  pulsar,  but  not  for  Vela.  It  has  been  hypothesized  that 
the  Vela  glitches  are  due  to  "corequakes"  in  an  hypothetical 
crystalline  neutron  core  (Pines  et  al.  1974).  Other  "glitch" 
hypotheses  are  discussed  by  Manchester  and  Taylor  (1977,  pp.  190-197). 


I 


The  electromagnetic  properties  of  neutron  stars  are  again  fixed 
by  the  properties  of  matter  throughout  the  star  and  its  response  to 
electromagnetic  fields.  These  fields  are  presumed  to  be  connected  with 
the  birth  of  the  neutron  star. 

Formation  of  Neutron  Stars 

It  is  generally  agreed  that  neutron  stars  are  formed  in  supernova 
explosions.  However,  both  observationally  and  theoretically  the  present 
evidence  to  back  up  this  hypothesis  is  small.  Only  in  the  case  of  the 
Crab  pulsar  is  the  association  with  the  Crab  Nebula  supernova  remnant 
(SNR)  unquestioned.  The  Vela  pulsar  is  almost  certainly  associated 
with  the  Vela  SNR,  but  only  in  a  few  other  cases  is  there  even  a 
questionable  association  of  pulsars  and  SNR's  (Manchester  and  Taylor 
1977,  p.  167).  A  recent  lunar  occultation  study  of  P0950+09  provided 
some  evidence  of  a  faint,  radio-emitting  SNR  surrounding  the  pulsar 
(Gopal-Krishna  1978).  There  are  several  factors  which  could  account  for 
the  failure  to  observe  a  significant  SNR-pulsar  correlation  even  if  it 
does  exist.  Pulsar  emission  is  believed  to  be  beamed  -  in  this  case 
about  80%  of  pulsars  would  not  be  observable  as  their  beam  would  not 
sweep  Earth  (assuming  a  random  beaming  direction  and  a  beam  width 
compatible  with  observation) .  This  beaming  and  the  impossibility  of 
detecting  very  faint  pulsars  could  explain  why  so  few  observed  SNR's 
have  had  pulsars  observed  in  them.  SNR's  generally  expand  and 
dissipate  with  age.  Thus  if  the  mean  pulsar  age  is  greater  than  one  to 
two  million  years,  any  SNR  which  may  have  surrounded  the  young  pulsar 
might  by  now  have  dissipated. 


The  estimated  occurrence  rate  of  supernovae  in  our  Galaxy, 
based  on  observation  and  on  stellar  evolution  theory,  is  between  one 
per  20  years  (based  on  SN  in  other  galaxies)  and  one  per  one  hundred 
years  (radio  SNR) .  This  is  in  accord  with  the  lowest  estimated  pulsar 
birthrates  (Manchester  1979).  It  has  been  suggested  (Cameron  1969) 
that  non-supernova  processes  may  also  form  neutron  stars. 

From  the  theoretical  point  of  view,  the  late  stages  of  evolution 
of  the  high-mass  stars  thought  to  be  supernova  progenitors  are  not  well 
understood.  Calculations  in  the  1960's  showed  that  only  extremely 
massive  stars  could  avoid  total  disruption  by  a  supernova  explosion  and 
so  leave  a  remnant.  This  suggested  that  neutron  stars  would  have  to  be 
formed  otherwise,  possibly  in  the  slow  contraction  of  a  slightly  over- 
massive  white  dwarf  (Cameron  1969) . 

More  recent  calculations  (Wheeler  1973;  Bruen  1973)  have 
indicated  that  for  certain  reasonable  predetonation  central  density 
regimes  (in  the  progenitor  star),  neutron  stars  should  form  in  supernova 
explosions  of  degenerate  stellar  cores  or  of  white  dwarfs. 

Massive  close  binary  systems  are  expected  to  evolve  very 
quickly  to  the  stage  where  one  of  the  stars  becomes  a  SN.  The 
material  of  the  other  star  accretes  onto  the  neutron  star  left  by  the 
explosion,  and  in  this  stage  the  system  is  an  X-ray  source  (possibly 
an  X-ray  pulsar).  As  this  star  evolves,  it  also  undergoes  a  SN 
explosion,  and  a  pair  of  high  velocity  neutron  stars  will  result,  or  a 
binary  radio  pulsar.  This  model  explains  observational  statistics  of 
the  various  guises  of  the  binary  system  (close  unevolved  binary, 
binary  X— ray  source,  pulsars)  quite  well  and  is  the  object  of  consider¬ 
able  study  at  present.  See  Manchester  and  Taylor  (1977,  p.  97)  for 


further  details  and  references. 


That  high  magnetic  fields  are  expected  in  neutron  stars  is  seen 
by  considering  that  they  are  the  result  of  the  collapse  of  stellar 
material  which  is  ionized  and  thus  a  good  conductor.  Thus  the  magnetic 
flux  (through,  say,  the  median  plane  of  the  star)  is  conserved,  so  that 

I  „  (^2  (2-4) 

where  BQ  and  B  are  the  magnetic  fields  at  corresponding  points  in  the 
star  before  and  after  collapse  (here  assumed  to  be  merely  a  scaling 
radially),  and  R0  and  R  respectively  the  total  radii  before  and  after 
collapse. 

A  contraction  factor  of  10 5 ,  which  is  that  needed  to  "shrink"  a 
star  to  the  size  of  a  neutron  star,  would  result  in  a  magnetic  field 
amplification  of  lO1^.  The  fields  of  500  to  3.4  x  104  G  observed  in 
magnetic  stars  could  be  "amplified"  to  the  magnitude  needed  (1013  G) 
for  current  pulsar  theories.  It  is  also  possible  that  magnetic  fields 
are  amplified  in  the  cores  of  massive  stars  before  a  SN  explosion 
(Angel  1978).  These  core  fields  would  not  be  observable  at  the  surface, 
and  thus  the  relatively  rare  magnetic  stars  would  not  be  the  only 
possible  neutron  star  progenitors.  The  fact  that  this  mechanism  would 
only  occur  in  the  most  massive  stars  would  explain  the  apparent 
paradox  (Ruderman  1972)  in  the  contraction  hypothesis  that  only  3%  of 
white  dwarfs  have  intense  magnetic  fields.  (A  relatively  small 
proportion  of  WD  progenitors  would  be  massive  enough.)  The  observa¬ 
tional  evidence  for  a  field  of  about  10 12  G  in  the  binary  X-ray  pulsar 


Her  X-l  has  been  cited. 

Differential  rotation  in  a  collapsing  star  could  cause  the 
generation  in  the  pulsar  interior  of  a  toroidal  magnetic  field  of  order 
10 1 5-10 1 7  G  (Pacini  1973).  The  stresses  produced  by  this  field  would 
cause  a  prolate  distortion  and  polar  wandering,  and  would  stabilize  the 
star  against  wobble  from  "glitches"  and  against  alignment  of  magnetic 
and  rotation  axes  if  they  are  oblique  (Goldreich,  Pacini,  and  Rees  1971). 
This  field  would  not  be  of  direct  significance  at  the  surface. 

A  similar  scaling  occurs  for  the  angular  velocity,  assuming  a 
simple  scaling  of  mass  (this  is  very  implausible  but  will  provide  a 
rough  estimate).  Amplification  of  angular  velocity  by  about  1010 
would  decrease  periods  from  105  to  106  sec  (main  sequence  stars' 
average)  to  10“5  to  10-1+  sec.  We  could  thus  expect  neutron  stars  to 
have  short  periods,  although  other  factors  (gas  ejection  in  SN,  the 
density  distribution  differences  of  progenitor  and  neutron  star, 
relativistic  constraints)  clearly  enter  a  more  quantitative  discussion 


of  this  point. 


CHAPTER  3 


NEUTRON  STAR  SURFACES 

The  fact  that  quantum-mechanical  effects  at  the  surface  of  a 
neutron  star  could  greatly  influence  magnetospheric  structure  was 
mentioned  in  the  classic  paper  of  Goldreich  and  Julian  (1969) .  Since 
then  these  effects  have  been  the  subject  of  controversy.  This  chapter 
summarizes  in  a  clear  and  consistent  manner  ideas  originally  put  forth 
in  a  series  of  confusing  papers  by  Ruderman,  alone  (1974,  1975),  with 
Sutherland  (1975),  and  with  Chen  and  Sutherland  (1974).  The  probable 
retention  of  ions  in  a  pulsar  surface,  as  predicted  by  the  (incorrect) 
calculations  of  Chen,  Ruderman,  and  Sutherland  (1974),  is  discussed. 

Effect  of  Magnetic  Fields 

In  the  interior  region  of  neutron  stars  the  effects  of  large 
magnetic  fields  are  negligible  compared  with  those  of  degeneracy 
(Canuto  1975b). 

However,  at  the  surface  the  fields  which  are  inferred  to  exist 
(of  order  1012-1013  G) ,  profoundly  affect  the  properties  of  matter,  in 
particular  the  opacity,  conductivity,  and  emissivity  of  ions.  The 
nature  of  the  surface  in  turn  determines  the  nature  of  the  magneto¬ 
sphere  (Ruderman  1975). 

The  opacity  due  to  the  three  major  contributors  -  the  free-free, 
bound-free,  and  Thomson  scattering  processes,  for  photons  with  E 
perpendicular  to  H  —  can  be  shown  (Lodenquai  et  al.  1974)  to  be  a 
factor  (w/a)H)2  smaller  in  the  presence  of  a  magnetic  intensity  H  than 
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magnetic  force  is  larger  than  the  nuclear  Coulomb  force  on  the 
electrons,  must  be  understood  in  order  to  specify  the  electrical 
properties . 

The  classical  relativistic  motion  of  a  charged  particle  in  a 
magnetic  field  of  strength  B  is  a  helical  orbit  whose  axis  parallels 
the  field,  and  whose  radius  is  the  Larmor  radius; 

R  _  IZIl!  .  (3-D 

L  eB  eB 

where  vj_  is  the  component  of  the  particle's  velocity  perpendicular  to 
the  field,  and  ymc2  the  particle's  energy.  RL  decreases  with  increasing 
B,  and  it  is  clear  that  when  is  of  the  same  order  as  the  particle's 
de  Broglie  wavelength,  the  motion  must  be  quantized.  In  similar 
fashion  to  the  quantization  procedure  for  the  Bohr  atom,  it  is  required 
that  the  action  in  the  plane  of  motion  be  quantized. 

The  old  quantum  theory  applied  to  a  spinless  electron  (Huang 
1963)  gives  the  quantization  condition  for  the  canonical  momentum  p 
(Landau  and  Lifshitz  1977,  p.  171) 

_  _  1 

§p*d£  =  (j  +  2)h,  j  =  0,1,...  (3-2) 

where  the  integral  is  around  one  orbit.  The  appropriate  Hamiltonian 
is 

H(p,r)  =  ^(P  +  f  A)  (3-3) 

where  e  is  the  magnitude  of  the  electron  charge.  The  equation  of 
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in  zero  field,  where  w  <<  am,  and  -from  =.  ti—  =  104-105  eV  for  a 

n  me 
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10  -10  G  surface  field.  An  heuristic  explanation  of  this  effect  for 

free-free  scattering  (Canuto  1975b)  is  that  the  electrons  cannot  inter¬ 
act  normally  with  the  E  field  when  it  is  forcing  them  to  move 
perpendicularly  to  magnetic  field  lines.  This  reduction  of  opacity 
speeds  the  cooling  of  neutron  stars.  However,  slowdown  of  the  internal 
superfluid  neutrons  is  a  source  of  thermal  dissipation  which  will  become 
important  at  a  later  stage  (age  106  yr)  of  a  pulsar’s  life  to  keep 
surface  temperatures  about  105-106  K  (Manchester  and  Taylor  1977, 
p.  175).  Young  pulsars  such  as  the  Crab  (103  a)  and  the  Vela  (104  a) 
could  have  temperatures  in  the  106-107  K  range  and  thus  emit  unpulsed 
X-rays  as  black-body  radiation  (Tsuruta  1974).  For  the  Crab, 
observations  failing  to  detect  such  X-rays  put  an  upper  limit  on  the 
surface  temperature  of  4.7  x  106  K  (Manchester  and  Taylor  1977,  p.  175). 
The  detection  of  unpulsed  optical  emission  from  both  the  Crab  and- Vela 
pulsars  (Peterson  et  al.  1978),  which  is  many  orders  of  magnitude  more 
intense  than  that  expected  from  a  neutron  star  emitting  as  a  black 
body  at  T  =  2  x  106  K  suggests  a  softening  of  higher-frequency 
radiation,  as  might  occur  in  the  case  of  X-rays  emitted  into  plasma 
surrounding  the  pulsars.  The  detection  of  unpulsed  X-rays  from  pulsars 
would  be  very  important  in  fixing  our  ideas  about  their  surface 
temperature. 

The  surface  electrical  properties  are  important  in  determining 
the  structure  of  the  neutron  star  magnetosphere  since  they  determine 
its  inner  boundary  conditions.  The  motions  of  the  surface  constituents 
in  intense  magnetic  fields,  and  the  properties  of  atoms  when  the 
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yields 


=  ®P. 


I  v  |  =  v  =  — —  B 
me 


(3-4) 


where  p  is  the  orbital  radius,  B  the  field,  and  m  the  electron  mass, 
v  being  the  electron  velocity. 

The  canonical  conjugate  of  the  p  in  equation  (3-2)  is  the 
position.  Its  derivative,  the  velocity  v,  is  thus  given  by  Hamilton's 
equation: 


v  =  VpH  =  ±<p  +  |  A) 


(3-5) 


so  that 


p  =  mv - A 

r  c 


(3-6) 


and  (3-2)  becomes 


§  (mv  -  A)  *d£  =  (j  +  2)h 
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(3-7) 
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The  first  part  integrates  directly  by  symmetry  to  2iTpmv  =  2irp^— ,  while 
the  second  has 


§A*d£  =  //V  x  A*ndS  =  // B*nds  =  iTp2B 


(3-8) 


giving 


TTp 


2  ^  =  (j  +  5)h 


or 


P2  -  ffU  +  2>* 
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(3-9) 


The  energy  is  the  numerical  value  of  H  in  such  a  system,  so 
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E  =  ta(P  +  cA)  =-|v 
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(3-10) 


To  describe  real  electrons  the  spin  energy  -y*B  can  now  be  added 

where 
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(3-11) 


and  spin  a  can  be  either  parallel  or  antiparallel  to  B.  Thus  the 
energy  levels  can  be  described  as 


E  =  —  nB  =  mc2n|  n=  0,1,2...  (3-12) 

me  Bq 

2  3 

where  Bq  (usually  denoted  Hq)  is  =  4.414  x  10 13  G. 

Including  a  possible  free  motion  along  the  magnetic  field,  with 
momentum  pz, 

o  B  P72 

E  =  mc2n^-  +  ^-  (3-13) 

Bq  2m 


is  the  non-relativistic  kinetic  energy  portion  (obviously  valid  only 
when  nB  <<  Bq)  of  the  exact  Dirac  equation  energy  set  found  by  Rabi 


(1928) 


E  =  me 2  [  1  +  ( — ) 
me 
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(3-14) 


The  level  spacing  in  the  non-relativistic  approximation  is 
mc2B/Bq  which  is  of  order  10  to  100  keV  in  the  fields  of  10 12  to  10 13  G 
of  a  canonical  magnetic  neutron  star.  It  is  likely  then  that  the  vast 
majority  of  electrons  are  in  the  magnetic  ground  state,  in  a  helix  of 
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radius  pQ  =  2  =  2.6  x  10-10/Bi2^  cm,  where  B12  is  in  units  of  1012  G. 

For  fields  of  order  10 12  G,  this  magnetic  (Landau)  radius  is  much  smaller 
than  the  Bohr  radius  aG  =  5  x  10~9  cm,  indicating  that  the  role  of  a 
proton's  Coulomb  field,  if  present,  would  be  that  of  a  perturbation. 

Since  nuclei  will  be  present  at  the  surface,  these  Coulomb  fields 
and  the  structure  of  atoms  in  very  strong  magnetic  fields  (B  >  10 10  G) 
should  be  considered.  A  light  atom  takes  the  shape  of  a  tube  of  the 
Landau  radius  with  a  length  determined  by  Coulomb  attraction.  Treatment 
by  perturbation  methods  (Landau  and  Lifshitz  1977,  p.  460)  or 
variationally  (Ruderman  1974)  gives  a  hydrogen  ground-state  energy 
(relative  to  motion  in  the  lowest  Landau  level  in  the  absence  of  a 
proton)  of 

E  =  fir? £n  =  -  2 a2mc2£n(^  (3_15) 

e2  1 

where  a  =  ^  is  the  fine  structure  constant.  The  quantity 

nc  Id  / 

1  2  2 

2a  mc  is  the  binding  energy  of  hydrogen  in  the  absence  of  a  magnetic 
field,  13.4  eV,  so  numerically 

E  ~  -13.4  in2 (19000^-) eV 

Bq 

^  -700  to  -1200  eV  for  B  =  1012-1013  G.  (3-16) 

These  values  are  reduced  somewhat  in  more  accurate  calculations,  but 
nevertheless  serve  to  show  that  very  large  binding  energies  result  from 
the  "squeezing"  of  the  atomic  radius  down  to  that  of  a  Landau  orbit, 
thus  causing  the  electron  to  move  deeper  into  the  protonic  Coulomb  well 
than  is  normally  possible.  The  length  of  the  atomic  "tube"  is  about 


a0/£n(a0/p0)  where  aQ  is  the  Bohr  radius. 


The  consideretion  of  multi-electron  atoms  is  facilitated  by 
solving  the  Schrodinger  equation  (Landau  and  Lifshitz  1977,  p.  459)  for 
wave  functions  associated  with  the  Landau  levels  of  a  spinless  electron, 
which  can  then  be  used  in  a  Hartree  calculation.  The  solutions  are 
(Cohen  et  al.  1970) 
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(3-17) 


where  to  (5)  is  the  normalized  confluent  hypergeometric  function 


F(-np, | m |  +  1,5)  (a  Laguerre  polynomial)  and  5 


=  „2  „  i  £ 


,  p  being 


2hc  2  pG 

the  radial  coordinate.  The  positive  integer  np ,  corresponding  to  the  n 
of  equations  (3-13)  and  (3-14),  denotes  the  Landau  level,  and  should  be 
zero  as  we  have  seen,  due  to  the  large  level  spacing.  The  new  quantum 
number  m  introduces  degeneracy  to  the  levels  because  fornp  =  0  the 
energy  eigenvalues  are  not  m-dependent.  The  motion  along  field  lines 
is  that  of  a  free  particle,  exp  [ipzz].  The  introduction  of  a  perturb¬ 
ing  Coulomb  field  will  remove  the  degeneracy  in  m,  and  also  quantize 
the  motion  along  the  field  lines.  Calculations  using  a  Hartree 
approximation  have  shown  that  the  filling  of  orbitals  in  multi-electron 
atoms  follows  m  rather  than  the  quantum  number  associated  with  the 
z-motion  (Chen  et  al.  1974).  The  Hartree  wave-function  for  an  atom  of 
atomic  number  Z  is  (Cohen  et  al.  1970) 
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(3-18) 


where  R^Pm)  =  Nmpmexp  (-  P 
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)  (since  n 


represents  the  z-dependence  of  the  m 


th 


p  =  0  for  all  m) ,  and  fm(zm) 
particle’s  wave  function.  The 


atom’s  energy  was  minimized  by  varying  fm(Zm)  through  a  parameter  am: 

u  .  . 

^m(^m)  =  (am)  exP  (“0^  |  Zm  |  ) 

In  this  way  the  total  binding  energies  of  atoms  and  their  once-ionized 
states  were  found.  The  ionization  energy,  which  then  follows  directly, 
is  of  order  135  eV  for  iron  at  2.2  x  1012  G. 

It  is  interesting  to  note  that  the  radius  where  an  electron  of 
np  =  0  and  given' m  is  most  likely  to  be  found  is 

pm  =  (2m  +  l)^Po  (3-19) 

and  is  thus  only  slightly  m-dependent.  The  electrons  of  a  many- 
electron  atom  represented  by  (3-18)  thus  form  a  "hollow  sheath"  around 
a  field  line. 

Such  cylindrical  atoms  bind  together  (Kadomtsev  and  Kudryatsev 
1971a)  end-to-end,  and  form  long  chains  consisting  essentially  of  a 
sheath  of  electrons  forming  a  Fermi  gas  in  the  "sheath,"  with  the  nuclei 
spaced  along  the  field  line  centred  within  the  sheath.  Such  chains  in 
turn  bind  to  one  another  to  form  a  lattice  when  alternate  ones  are 
displaced  by  one-half  the  inter -nuclear  spacing  (Chen  et  al.  1974). 

Emission  Properties 

For  the  mainly  56Fe  composition  expected  to  be  characteristic 
of  the  surface  of  neutron  stars,  in  the  magnetic  field  of  2  x  1012  G, 
the  ions  will  have  a  cohesion  energy  of  4.1  keV  (Flowers  et  al.  1977), 
and  the  electrons  in  the  sheath  will  have  a  750  eV  Fermi  energy  (Chen 
et  al.  1974).  The  incorrect  prediction  by  Chen,  Ruderman,  and 
Sutherland  (1974)  that  the  cohesive  energy  of  atoms  to  chains  would  be 
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14  keV  at  2  x  10 12  G  led  them  to  conclude  that  this  binding  into  chains 
would  alone  prevent  emission  of  ions  from  the  surface.  The  lower 
cohesive  forces  now  calculated,  while  implying  that  cohesion  may  not  be 
entirely  responsible  for  ion  retention,  do  not  imply  that  ions  flow 
freely  from  the  surface.  The  atoms  which  may  be  removed  thermally  from 
chains  will  form  but  a  small  percentage  of  the  total  number  of  atoms  at 
the  pulsar  surface,  since  kT  -  10  eV  <<  cohesive  energy.  Thus  the 
Fermi  energy  of  750  eV,  mentioned  previously,  will  remain  unchanged  or 
nearly  so.  In  any  ionization,  itself  requiring  135  eV,  the  electron 
would  have  to  acquire  more  than  750  eV  to  fill  an  allowed  level  above 
the  Fermi  sea.  Due  to  this,  free  atoms  on  a  pulsar  surface  cannot  be 
ionized  except  by  extremely  energetic  X-rays,  which  would  probably  have 
a  very  low  interaction  cross  section  with  the  outermost  electrons  of 
any  atom  present.  This  result  would  apply  not  only  to  iron  but  also  to 
other  atoms  (such  as  helium  (Michel  1975b))  which  may  be  present  in 
trace  amounts  (see  Fawley  et  al.  1977).  The  fact  that  ionic  field 
emission  would  be  suppressed  due  to  the  fact  that  neutron  star  surfaces 
are  protected  against  ion  formation  (from  atoms) ,  due  to  their  high 
Fermi  energy,  has  not  been  discussed  in  the  literature. 

It  is  worth  noting  that  this  protection  mechanism  is  similar  to 
that  which  allows  neutrons  to  exist  stably  in  the  interior  of  a  neutron 
star.  Ruderman  and  Sutherland  (1975)  estimate  that  an  applied  field 
E  such  that  Zce£E  is  the  order  of  the  cohesive  energy  could  separate 
ions  of  charge  Zce,  consisting  of  iron  nuclei  and  several  inner 
electrons.  Even  with  the  reduced  cohesive  energies  recently  calculated, 
a  field  of  about  10 11  V/cm  would  be  required.  While  a  field  of  this 


order  could  be  obtained  for  a  pulsar  rotating  in  vacuum  (see  Chapter  4) , 
screening  will  occur  in  a  realistic  pulsar  magnetosphere  to  prevent  such 
high  fields  from  developing.  Ion  retention  thus  will  characterize 
pulsar  surfaces. 

Ruderman  and  Sutherland  (1975)  further  argue  that  thermionic 
emission  and  sputtering  due  to  the  impact  of  relativistic  electrons 
and  positrons  are  negligible  (however,  see  Chapter  4). 

The  emission  of  electrons  depends  very  sensitively  on  the  work 
function  for  the  surface.  This  is  not  well  determined  but  should  be  of 
the  order  of  one  keV.  However  electron  emission  is  probably  much  easier 
than  that  of  ions  (Ruderman  1972a). 

The  conductivity  of  the  atomic  chains  is  essentially  infinite 
parallel  to  the  chain  axes  (field  lines),  since  the  electrons  form  a 
one -dimensional  Fermi  gas  (Ruderman  1974).  In  the  perpendicular 
direction  it  will  be  less  (Ruderman  1971).  Thus  the  surface  conductiv¬ 
ity  is  anisotropic. 

Due  to  the  formation  of  atomic  chains  bound  to  one  another,  the 
surface  itself  will  end  abruptly,  and  be  more  similar  to  that  of  a 
planet  (or  of  "a  not  very  smoothly  shaved  porcupine,  with  chains  of 
Fe-atoms  sticking  out"  (Borner  1973))  than  that  of  a  star. 

The  gravitational  forces  at  the  surface  are  very  large,  and  for 
a  neutron  star  of  M  solar  masses  and  radius  Rg  (in  units  of  !06  cm)  the 
gravitational  potential  energies  of  electrons  and  protons  are  -8  x  104 
M/Rg  eV  and  -1.4  x  108  M/Rg  eV,  respectively  (Goldreich  and  Julian  1969). 
For  M  -  1  and  a  typical  Rg  =  1  these  are  about  -.1  MeV  and  -100  MeV 
respectively. 


To  summarize,  the  surface  of  a  rapidly  rotating  neutron  star  is 
subjected  to  the  action  of  a  very  strong  electric  field  sufficient  to 
remove  either  electrons  or  ions .  This  particle  outflow  will  guarantee 
that  the  star  has  a  magnetosphere  of  charged  particles,  and  the 
difference  in  ease  of  removal  of  electrons  and  ions  has  entered  recent 
models  in  a  very  significant  role  (Ruderman  1975). 


CHAPTER  4 


PULSAR  ELECTRODYNAMICS 

The  problem  of  the  pulsar  magnetosphere  has  not  yet  been 
"solved,"  in  the  sense  that  the  behaviour  of  fields  and  particles  in  the 
vicinity  of  a  rotating  magnetic  neutron  star  have  not  been  detailed  in  a 
self-consistent  manner.  Most  theories  of  pulse  production  consider  the 
magnetic  symmetry  axis  and  that  of  rotation  to  be  non-aligned  -  this 
departure  from  symmetry  presumably  being  reflected  in  the  highly  aniso¬ 
tropic  emission  pattern  thought  to  cause  the  pulses.  However,  such  a 
problem  is  considered  intractable,  and  most  effort  has  gone  into  the 
simplified  case  of  parallel  or  antiparallel  magnetic  and  rotational 
axes.  This  special  case  removes  both  implicit  time -dependence  and  axial 
asymmetry,  and  presumably  the  pulse,  from  pulsars.  Despite  these 
simplifications,  even  this  problem  has  not  been  satisfactorily  solved, 
partly  because  the  magnetospheric  plasma  still  retains  many  degrees  of 
freedom,  and  also  due  to  an  incomplete  understanding  of  the  extreme 
physical  conditions  and  of  the  boundary  conditions  involved. 

The  aligned  rotator  (magnetic  and  rotation  axes  parallel  or  anti¬ 
parallel)  will  be  considered  first  in  some  detail.  The  general  features 
it  possesses  should  be  understood  before  we  very  briefly  consider  the 
more  general  non-aligned  (oblique)  rotator,  to  which  they  may,  to  an 
extent,  also  apply.  Numerical  models  of  non-aligned  rotators  are  of 
some  interest.  After  the  basic  structure  of  the  magnetosphere  is 
explained,  in  this  chapter,  the  emission  mechanisms  which  may  account 
for  the  pulsed  radiation  will  be  considered,  in  the  next. 
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In  this  chapter,  unless  otherwise  specified,  the  neutron  star 
itself  will  be  considered  as  a  rigidly  rotating  perfect  conductor. 

While  this  approximation  breaks  down  in  some  regions  of  the  star,  it 
will  represent  the  overall  steady-state  situation  accurately. 

That  electric  fields  play  a  major  role  in  determining  the  magneto- 
speric  structure  was  first  suggested  by  Goldreich  and  Julian  (1969) ,  for 
the  aligned  rotator.  The  same  basic  argument  holds  for  most  of  the 
surface  of  an  oblique  rotator  (Cohen  and  Toton  1971) . 

The  Vacuum  Rotator 

A  perfect  conductor,  rotating  in  a  magnetic  field  in  a  vacuum, 
acquires  a  charge  distribution  in  its  interior,  which  in  turn  produces 
an  electric  field.  In  a  pulsar  the  force  exerted  by  this  field  at  the 
surface  will  remove  charged  particles  to  form  a  magnetosphere.  The 
details  of  the  vacuum  aligned  rotator  problem,  which  follow,  are  as  in 
Israel  (1974).  Since  the  conductivity  is  perfect,  no  charged  particle 
within  the  star  can  feel  a  Lorentz  force.  At  all  internal  points 

E  +  -  x  B  =  0  (4-1) 

.  c 

The  particles  rotate  in  the  star  and  thus  have  velocity  v  =  w  *  r. 

Fixing  ui  =  wz  and  setting  up  non-rotating  (see  Webster  and  Whitten  1973) 

spherical  coordinates  centred  on  the  star,  and  with  the  magnetic  field 

9  9 

symmetric  about  u>,  the  rotational  symmetry  means  that  =  0. 

_  3^ 

Immediately  follows  since  E  =  -V$  -  — ,  where  $  is  the  scalar,  and 
A  the  vector  potential.  The  other  components  of  E  are  most  easily 

—  A 

expressed  by  (4-1)  where,  now,  v  =  mrsin6^; 


(4-2) 


Eg  =  -(t-  x  B) q  =  -^torsinQ  Br 
c  c 


Er  =  -(—  *  B)r  =  -^mrsin6  Bg 


Once  the  magnetic  field  components  are  specified,  these  relations,  true 
within  the  star,  allow  the  induced  charge  distribution  to  be  calculated 
using 


p  =  —  VE  = 


z-(B  "  T  r  x  (V  x  B)} 
2ttc  2 


(4-3) 


This  induced  charge  distribution  gives  rise  to  an  electric  field 
externally  as  well  as  internally. 

The  external  fields  may  be  found  by  noting  that  in  the  vacuum 

surrounding  the  star  Maxwell's  equations  are  homogeneous  and  thus  both 

the  electric  and  magnetic  fields  are  derivable  from  potentials.  These 

potentials  will  be  fixed  by  boundary  conditions.  In  this  case  the 

tangential  component  Eg  of  E,  and  the  normal  component  Br  of  B  are 

continuous  at  the  surface  (r  =  a) .  The  electrostatic  potential  satisfies 

$  -  0(-^r)  as  r  -*  00 ,  and  solves  V2$  =  0. 
rz 

In  the  case  of  a  dipole  magnetic  induction  field  of  moment  y, 


Br  =  -^jcos0  =  B0(-p)  cos6,  and  B0  =  ™9-  =  ~y(r^)  sin9 


(4-4) 


where  BQ  =  Br  (r  =  a,  6  =  0)  =  is  the  magnetic  induction  (loosely, 

the  magnetic  field)  at  the  pole.  The  boundary  condition  on  the 

continuity  of  Eg,  from  (4-2)  and  (4-4),  is  Eg  =  -  ~('|f')r=a 
B  3  * 

- - 2.  M(a  cosQsing. )  which  restricts  the  solution  of  Poisson's 

r  r^  r-d 


equation  for  the  electrostatic  potential  to 


$(r  >  a,0)  =  - 


l  B0wa5  P2(cos0) 


(4-5) 


3  c  r 3 


The  field  in  (4-4)  does  not  include  some  source  of  field  (a 
point  dipole),  which  would  appear  as  a  singularity  at  the  origin  (see 
e.g.  Jackson  1975,  pp.  138-141).  Mestel  (1971)  notes  that  (4-4)  would 
imply  that  the  curl  of  the  magnetic  induction  would  be  everywhere  zero, 
which  is  unrealistic.  Within  the  star  he  suggests  that  the  field  be 
regarded  as  uniformly  B0z.  The  discontinuity  in  Bg  outside  and  inside 
would  produce  a  surface  current  (Jackson  1975,  pp.  19-21),  in  the 

A  A 

direction  f  x  0  =  <J>.  Mestel  s  model  also  gives  a  different  surface 
charge  distribution. 


(4-6) 


a 


than  will  be  found  in  (4-8)  from  Goldreich  and  Julian.  In  light  of  the 


anisotropy  of  the  surface  conductivity,  as  discussed  in  Chapter  3,  it 
is  worth  noting  that  the  Goldr eich-Julian  solution  depends  only  on  the 
Eg  boundary  condition.  With  infinite  conductivity  in  only  the  0 
direction  this  boundary  condition  will  not  change  despite  the  fact  that 
the  conductivity  in  r  and  <J>  may  be  very  small. 


However,  in  this  case  the  continued  existence  of  a  current  in 


the  toroidal  direction  would  be  questionable.  This  problem,  and  that 
of  an  anisotropic  oblique  dipole  rotator  (discussed  later),  would  be  of 


great  interest. 

The  external  electric  field  from  (4-5)  has  a  radial  component 


Er  -  -!i 

r  3r 


BO0)ab  1 


cr 


—  2-(3cosz0  -  1) 


BQu)a 

2c 


(3cos20  -  1)  (r  =  a) 


(4-7) 


while  the  radial  component  of  the  internal  field  is 


Bpcja3 

cr2 


-j  sin20 


Bpcja 

2c 


sin20 


(r  =  a) 


(4-8) 


The  induced  surface  charge  a  is  —  times  the  difference  of  (4-7)  and 

4tt 

(4-8),  that  is 


CT  =  "  cos2e  (4"9) 

More  important,  however,  is  the  fact  that  there  is  a  component  of  the 
electric  field  which  can  remove  charged  particles  from  the  surface.  In 
the  strong  magnetic  fields  near  the  surface  the  motion  of  charged 
particles  is  free  only  along  the  magnetic  field  lines.  Thus  the 

_  A 

component  of  E  along  B,  of  magnitude  E*B,  is  the  source  of  energy 
increases  of  charged  particles.  (E  *  B  drift  velocities  will  be  small 
since  |e|<<|b|  (Jackson  1975,  p.  582),  and  are  neglected.)  E*B,  from 
(4-4)  and  the  gradient  of  (4-5),  is 

2.  8 

E*B  =  -  -  cos30  (4-10a) 

cr  7 

Inside,  that  E*B  =  0,  which  is  implied  by  the  Lorentz  Force  being  zero, 
may  be  easily  shown. 

The  problem  of  a  perfectly  conducting  oblique  rotator  in  vacuum 
was  solved  by  Deutsch  (1955).  For  a  dipole  field  (Jackson  1978b; 


' 


(correction  of  his  eq.  2.11)), 


E«B 


Bn2u)a8  1  r  ,rs  2 


cr 


—  4[(a}  (cos 


+  4  cos0cos20 


m 


cos0m  -  cos0  -  4qs(^a3B°)) 
-  COS  X  COS0m  +  COS0] 


(4-10b) 


where  x  is  the  obliquity,  0  the  polar  coordinate  from  the  rotation 

axis,  and  0m  the  (time-dependent)  polar  coordinate  from  the  magnetic 

axis.  This  result  is  valid  only  near  the  star  (i.e.,  r  <<  — ) ,  and  a 

co 

possible  stellar  charge  qs  has  been  included.  If  qs  =  0,  the  difference 
from  (4-10a)  lies  only  in  the  time-dependent  angular  factor  of 
magnitude  about  unity.  Finding  the  effect  of  the  anisotropic 
conductivity  at  the  surface  of  a  neutron  star  on  Deutsch's  vacuum 
solution  would  be  an  interesting  problem  in  electrodynamics.  Deutsch's 
general  solution  has  surface  currents  which  flow  in  both  the  poloidal 
(0)  and  toroidal  (cf>)  directions.  At  the  neutron  star  surface  only 
poloidal  flow  would  be  resistance-free.  There  would  undoubtedly  still 
be  a  similar  value  of  E*B,  although  the  angular  dependence  would  be 
different  from  that  of  equations  (4-10). 

The  important  result  is  that,  depending  on  the  surface  condi¬ 
tions,  positive  or  negative  charged  particles  may  be  removed  from  a 
pulsar  surface.  A  rotating,  magnetized  neutron  star  will  not  be 
surrounded  by  vacuum  (Goldreich  and  Julian  1969),  but  rather  by  a 
magnetosphere,  or  magneto-electric  atmosphere  (Jackson  1978a),  formed 
of  charged  particles  extracted  from  the  star. 

The  rate  at  which  particles  are  pulled  from  the  surface  will 
obviously  depend  on  their  binding  energy  or  work  function.  As  discussed 
in  Chapter  3  this  differs  for  electrons  and  for  ions  and  may  be  of  large 


magnitude  (Ruderman  1975) . 

Another  mechanism  for  "populating"  the  vicinity  of  the  neutron 
star  is  the  production  of  electron-positron  pairs  by  Y-raYs  moving  in 
the  intense  magnetic  fields  (Sturrock  1971;  Erber  1966) .  The  creation 
of  pairs  from  the  energy  of  the  magnetic  field  alone  is  not  possible 
(Chiu  1969). 

However,  before  investigating  the  effects  of  a  neutron  star's 
plasma  atmosphere,  the  radiation  from  an  oblique  magnetic  dipole 
rotating  in  vacuum  should  be  considered,  as  a  rough  guide  to  possible 
radiation  mechanisms. 

The  power  radiated  by  a  non-constant  magnetic  dipole  y  is  given 
by  (Landau  and  Lifshitz  1962,  p.  217): 

dE  =  _  2 
dt  3  c3 


Here  y  =  io2ysinx,  where  the  dipole  makes  an  angle  x  with  the  angular 

-  g  a3 

velocity  5.  The  rotating  frame  value  of  y,  — ^ — ,  will  be  about  the 

same  in  the  non-rotating  system  which  we  are  using  due  to  the  fact 

that  aco  <<  c,  so  that  the  power  radiated  will  be 


dE 

dt 


1  q)4a6Bo2sin2x 


(4-11) 


dE 

However,  for  a  rotating  body  losing  energy,  —  =  Iwm,  where  I  is  the 
moment  of  inertia. 


Equating  these  expressions  for  energy  loss  will  give  the  expected 


rate  of  change  of  rotation  frequency  for  a  magnetic  dipole  associated 
with  a  spinning  body: 


(4-13) 


-1  .  -2 . 

or,  since  P  =*  2irw  and  P  =  -2tto)  w, 

p  ,.l  a6BQ2sin2X  p-i 
24tT2  Ic3 

That  all  period  derivatives  are  positive  is  heartening  at  this  point 
(it  seems  to  at  least  confirm  that  a  spinning  body  is  involved) ; 
however  the  fact  that  there  is  no  marked  correlation  of  P  and  £  implies 
that  if  the  basic  power-loss  mechanism  is  magnetic-dipole  radiation, 
the  parameters  a,  B0,  siny,  and  I  vary  widely  from  pulsar  to  pulsar. 

For  a  "canonical"  pulsar  with  siny  =  1,  a  =  108  cm,  and  I  =  1045  gem2, 
we  have  (Manchester  and  Taylor  1977,  p.  177),  for  P  in  seconds, 

B0  =  3.2  x  1019(PP)^  G  (4-14) 

which  gives  values  between  about  2  x  10 10  G  and  2  x  10 13  G,  with  10 12  G 
being  typical.  Using  these  "typical"  values  in  (4-11),  one  finds  that 
by  this  mechanism  the  Crab  pulsar  would  put  out  about  7  x  1038  erg-s-1, 
which  is  just  about  what  is  needed. 

A  further  insight  into  pulsar  emission  processes  can  be  gained 
by  hypothesizing  that  for  a  pulsar  with  given  parameters  (which  fix  the 
constant  K) ,  the  torque  slowing  rotation  will  be  proportional  to  a  power 
of  the  rotation  frequency,  so  that 

(L  =  -Km11  (4-15) 

n  being  called  the  braking  index.  For  magnetic  dipole  radiation  we 
have  seen  that  n  =  3.  Differentiation  of  (4-15)  and  elimination  of  K 


give 


n  = 


0)0) 


(4-16) 


so  that  the  braking  index  can  be  determined  from  timing  results  if  they 
are  good  enough  to  give  a  second  period  (or  frequency)  derivative.  For 
example  the  braking  index  of  the  Crab  Pulsar  P0531+21,  calculated  from 
recent  timing  results  (Groth  1975),  is  2.5071.  This  has  been  seen  as 
indicating  that  the  magnetic  dipole  model  serves  as  a  good  starting 
place  in  constructing  pulsar  models.  It  is  rather  disturbing  that  second 
derivatives  recently  published  for  9  other  pulsars  (Gullahorn  and  Rankin 
1978b)  lead  to  very  large  braking  indices  (Appendix  2) . 

If  the  constant  K  in  equation  (4-15)  is  indeed  constant  through¬ 
out  the  life  of  a  pulsar,  (4-15)  may  be  integrated  to  give  an  expression 
for  the  pulsar  age  t: 


t 


0) 


(n-l)d) 


03  n_1 

[1  -  (— )  ] 

03  n 


(4-17) 


where  03o  is  the  ’’birth"  rotation  frequency.  If  03  <<  03o  this  age  is 
approximately  equal  to  the  characteristic  age  for  index  n: 


03  =  P 

(n-l)d)  (n-l)P 


(4-18) 


Our  previous  definition  of  characteristic  age  as  PP  1  is  simply  x 
for  n  =  3. 

A  simple  argument  (Canuto  1975c)  sbout  the  magnitude  of  the 
"birth”  rotation  frequency  is  that  it  is,  at  maximum,  that  frequency 
corresponding  to  breakup  (against  Newtonian  gravity)  by  centripetal 


forces,  so 


GM  v2  o 
— 7  =  —  =  u)0  a 
az  a  ° 


or 


w02  =  Gp 


(4-19) 


where  p  -  Ma  3 .  At  a  density  of  10 15  g  cm  3  and  for  a  =10  km,  to0  would 
be  104  rad  s-1,  so  that  even  for  the  Crab  (w  -  33  x  2tt  rad  s-1),  u)<<0)o. 


Canuto  points  out  that  expansion  velocities  observed  in  the  Crab  and 
Vela  supernova  remnants  would  also  be  consistent  with  an  w0  of  this 
size. 


The  above  considerations  are  for  a  rotator  in  vacuo,  but  the 


creation,  due  to  the  parallel  component  of  the  electric  field,  of  a 
plasma  atmosphere  about  the  neutron  star,  will  change  the  equations  to 
be  satisfied  in  the  vicinity  of  the  star.  In  particular  source 
densities  in  Maxwell’s  equations  will  be  non-zero.  The  solution  to  the 
entire  problem  will  then  differ  from  what  has  been  described.  The 
magnitude  of  the  source  densities  will  determine  how  much  the  models 
differ  from  the  vacuum  oblique  rotator. 

Magnetospheric  Equations 

To  consider  the  effects  of  source  densities,  the  following 
discussion,  based  on  Cowling  (1976)  and  Israel  (1974)  is  offered. 

The  field  equations  with  sources  are 


(4-20) 


V  x  E 


1  _3B 
c  3t 


(4-21) 


V*B  =  0,  V*E  =  4ttp 


(4-22) 


and  the  generalized  Ohm's  Law  is 


j  =  a(E  +  ^  x  b)  (4-23) 

Equating  the  currents  in  (4-20)  and  (4-23),  taking  the  curl,  and  using 
(4-21)  and  (4-22)  to  simplify  the  resulting  V  x  (v  x  B)  and  V  x  E 
terms,  one  obtains 

8B  —  —  i  2 

—  =  V  X  (v  X  B)  +  n(V2  -  ^2  B  (4-24) 

c2 

where  q  =  is  the  magnetic  diffusivity.  Ignoring  the  second-order 

time  derivative  of  B,  which  is  due  to  displacement  currents,  should  be 
valid  where  32  <<  1  (Jackson  1975,  p.  471).  In  any  case,  the  plasma 
around  a  neutron  star  would  be  expected  to  have  a  very  high  conductivity, 
so  that  q  <<  1  and  the  second  term  may  be  neglected.  The  resulting 

H  =  V  X  (V  X  B)  (4-25) 

may  be  shown,  by  considering  the  flux  traversing  any  moving  surface 
in  the  plasma  (Alfven  and  Falthammar  1963,  p.  101),  to  imply  that  the 
magnetic  field  is  "frozen  into"  the  plasma.  By  using  (4-21),  the  fact 
that  this  implies  that  the  Lorentz  force  is  zero  (the  "hydromagnetic 
condition"),  may  be  easily  noted.  Obviously  a  way  of  checking  models 
based  on  a  vanishing  Lorentz  force  ("force-free  magnetospheres")  would 
be  to  check  that  the  conditions  in  the  model,  once  solved,  are  such  that 
the  second  term  in  (4-24)  is  negligible. 

The  case  of  near  perfect  conductivity  and  a  magnetosphere 
corotating  with  the  neutron  star  due  to  the  freezing  in  of  the  magnetic 


field  will  now  be  considered. 
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The  structure  of  the  magnetosphere  near  the  neutron  star  will 
now  be  considered,  based  on  equations  given  in  Scharlemann  and  Wagoner 
(1973)  and  Brooks  (1974).  A  steady-state  solution  for  a  time- 
independent  aligned  rotator  will  be  sought,  with  the  assumptions  that 
only  electrons  and  a  single  ion  species  are  present  with  well-defined 
bulk  velocities,  and  that  gravitational  and  inertial  forces  may  be 
neglected.  Complications  which  make  the  structure  of  the  outer 
magnetosphere  non-trivial  are  mentioned,  but  only  the  near  magnetosphere 
is  considered  in  detail. 

In  the  inertial  frame,  the  Lorentz  force  has  form: 


-en_[E  +  (3-  x  B)  =  m_c2n_3_* V (y3) _  (4-26) 

Zen+[E  +  (3+  *  B)  =  m^c  2n+3+- V  (y3)  +  (4-27) 


so  these  are  the  equations  of  motion  of  electrons  and  ions  of  number 

densities  n_  and  n_j_,  masses  m_  and  m_|_,  and  velocities,  v+  =  3+c 

g 

respectively.  The  derivative  -r—  is  zero,  so  the  RHS  is  the  "inertial 

d  t 

force  term"  of  the  convective  derivative  m^-  =  m(-^-  +  v*Vv). 

dt  dt 

Neglecting  the  inertial  force  term  is  equivalent  to  assuming 

B2  o 

massless  particles.  This  may  be  realistic  where  —  <<  yyc^  (y  =  mass 
density) ,  and  in  turn  implies  infinite  conductivity  (collisions  are 
neglected),  so  that  either  algebraically,  or  arguing  that  the  Lorentz 
force  must  vanish  for  the  case  of  infinite  conductivity,  we  obtain  the 
condition  mentioned  after  (4-25), 


E  =  -3+  x  B 


(4-28) 


Neglect  of  inertial  force  terms  also  brings  in  an  indeterminacy 
of  electron  and  ion  velocities  -  only  the  current  j  =  e (Zn+8+-n_8_) , 
which  is  a  linear  combination  of  the  desired  velocities,  is  specified 
(Scharlemann  1974) . 

Maxwell's  equations 

V  •  B  =  0 - ►  B  =  V  x  A  (4-29) 

V*E  =  4rp  =  4ire  [Zn+-n_]  (4-30) 


are  independent  of  assumptions,  while 


V  x  B  =  -§  +  —  j 
c  9t  c 


in  this  case  becomes 


(4-31) 


V  x  B  =  4ire  (Zn+6+-n_3_) .  (4-32) 

The  remaining  Maxwell  equation: 

V  x  E  =  -  (4-33) 

implies  that  V  x  E  =  0,  so  that  E  =  -V$ . 

Equation  (4-28),  as  we  have  seen,  implies  that  E*B  =  0:  the 
magnetic  field  lines  are  equipotentials  if  (4-28)  holds.  It  has  already 
been  shown  to  imply  that  the  magnetic  field  lines  are  frozen  into  the 
plasma . 

Equation  (4-32),  crossed  into  B  gives 


(V  x  B)  x  B  =  471-6(211+8+  -  n_(3_)  x  B 


(4-34) 


which  using  (4-26)  and  (4-27)  becomes 


(V  x  B)  x  b  =  -47re(Zn+  -  n_)E  =  -(V*E)E 


(4-35) 


so 


°r  1£  it  *  0 


(V*E)E+(V  x  b)  x  B  =  0 


(V-E)E  +  (V  x  B  -  -  ^)  x  B  =  0 

C  dt 


(4-36) 

(4-36a) 


(Mestel  1973).  This  is  the  equation  of  a  force-free  electromagnetic 


field. 


The  axially  symmetric  corotating  magnetosphere  (Goldreich  and 


Julian  1969)  has  no  toroidal  electric  field  component  for  the  same  reason 

3  9 

that  the  aligned  vacuum  dipole  did  not,  assuming  —  =  —  =0.  (A  looser 

d(p  dt 

9  9 

condition,  —  =  -orr- ,  is  considered  by  Mestel  (1971,  1973).) 

O  t  o  (p 

E  =  -(3  x  B)  from  (4-28),  and  if  =  0,  (4-25)  shows  that  8  x  B  is 
<p  <f)  dt 

curl-free  (solenoidal) .  Writing  the  vectors  as  the  sum  of  poloidal 

_  «— •  A  A 

(8p,Bp)  and  toroidal  (E^^jB^)  vectors,  one  obtains  (Mestel  1961) 


V  x  (gp  x  Bp)  =  0  (4-37) 

and 

v  x  (e^i  x  Bp  +  8p  x  B^)  =  0.  (4-38) 

3p  x  Bp  is  a  toroidal  vector  (it  is  actually  E^),  and  being  solenoidal, 
must  be  the  (^-derivative  of  a  scalar  function,  which  by  axial  symmetry 
is  zero.  Thus  8p  x  Bp  =  0,  or  8p  =  KBp  where  K  is  a  scalar.  Since 
Bp  =  B  -  B^cj),  8  may  be  written  8  =  K(B  -  B^cj))  +  a$,  which  may  be  put 
into  (4-28)  to  yield 


E  =  (KB^  -  a$)  x  Bp. 


(4-39) 


corsinO 


To  determine  the  scalar  a  we  note  that  3  =  3^$  +  3p 
,  n  cjrsinQ  _ 

where  3^ - - -  for  corotation  at  radius  r.  Then  a  =  K3^  + 

so  that 

3  =  KB  +  -r-s-ln9-  $  =  KB  +  ~  X  r  . 

c  c 


(4-40) 


We  are  now  in  a  position  to  find  E  and  the  charge  density  which  it 
implies.  Since  E  =  -3  x  B,  the  divergence  of  E  is 


V*E  =  -V  (— -*  -r  x  B)  =  -B-V  x  (^  X-r)  +  ^r-s..ln9^.  (v  x  B)  (4-41) 

C  CO 

In  a  non-charge-separated  atmosphere  where  3+  =  3_  in  bulk,  or  in  a 
highly  charge-separated  region  where  either  n+  or  n_  may  be  neglected, 
(4-32)  and  (4-30)  imply  that 


V  x  b  =  3V*E  (4-42) 

and  expanding  V  x  (w  *  r),  noting  w  =  rnz, 

(V-E)  [1  -  -r-S-in-6^3]  =  -  |  [w(V-r)  -  (u*r)r]  =  -2B-w/c  (4-43) 

so  that 

V •  E  B’w  r,  u)2r2sin2e-,-1  ,,  //N 

p  =  TT  "  -  2^  [1 - ^ - ]  (4-A4> 

This  is  the  charge  density  derived  by  Goldreich  and  Julian  for  a 
perfectly  conducting  corotating  magnetosphere.  Corotation  cannot  hold 
out  to  the  light  cylinder  (the  locus  where  wrsin©  =  c) ,  and  indeed  the 
toroidal  current  associated  with  the  space  charge  in  (4-43)  becomes  a 
major  source  of  field  near  the  light  cylinder,  distorting  the  dipole 


field. 


Neglecting  this  distortion  one  can  delineate  the  zone  of 


corotation  roughly  by  that  dipole  field  line  which  first  touches  the 
light  cylinder.  A  line  touching  the  stellar  surface  at  colatitude  0Q 
has  a  parametric  equation  (Lorrain  and  Corson  1970,  p.  69;  Jackson  1976) 


r 


.sin  0 
a 

sxn0o 


(4-45) 


Setting  r  =  c/w  at  sin©  =  1  gives  the  value  sin0o  =  C^y-)^  and  allows 
the  polar  equation  of  the  edge  of  the  CRM  to  be  given; 


r  CRM  " 


csin20 

0) 


(4-46) 


The  surface  region  where  0  <  0Q  defines  a  "polar  cap"  of  the  pulsar, 
of  radius  (Ruderman  and  Sutherland  1975) 


rD  =  a(— )**  *  1.4  x  104  P-1  cm. 
P  v  c  ' 


(4-47) 


Mestel  (1971)  has  pointed  out  the  conditions  to  be  found  in  such 
a  corotating  magnetosphere.  The  charge  density  is  the  excess  of 
positive  over  negative  charges,  but  the  actual  situation  is  believed  to 
be  one  of  large  charge  separation.  This  is  equivalent  to  the  assumption 
that  electric  fields  due  to  space  charge  densities  are  negligible  in 
comparison  to  the  field  implied  by  the  force-free  condition.  In  this 
case  the  ratio  of  the  (excess)  charge  density  to  that  of  positive 
charge,  |p|/n+Ze,  will  be  at  least  of  order  unity.  In  contrast,  this 


ratio  in  a  normal  plasma  would  be  far  less  than  unity. 

Bu) 


Using  (4-44)  in  the  form  |p 
Pm 


2ttc 


,  and  noting  that  the  ionic 


number  density  n+  =  ,  where  A  is  the  mass  number  of  the  ions,  m^  the 

mass  of  a  hydrogen  atom,  and  pm  the  mass  density, 


0) 


0) 


-lei  K  r.PrnZe,  1 

n+Ze  2irc  L  AmjjJ 


(  B2  \  /  a)  \  _  ,  B2 
''2iTpinc2  ^ZeB/Am{jc' 


87rpmc2  ui 


(4-48) 

where  cjo-j_  =  is  the  gyration  frequency  of  the  ions.  Since  will  be 

of  order  10-16,  one  is  assured  that  the  macroscopic  equations  (4-30)  and 
(4-32)  are  valid.  We  have  already  concluded  that  the  density  of  the 
plasma  is  such  that  the  left-hand  side  of  (4-48)  is  of  order  unity. 

Then  the  ratio  of  magnetic  to  rest  mass  energy, 


B2 

8upc2 


0  (— “•)  >  >  1 . 

0) 


(4-49) 


As  seen  in  Chapter  3  this  implies  that  the  particle  motions  are 
relativistic.  Mestel  also  shows  that  the  convection  current  pv  is 
effectively  the  total  current,  and  that  the  electrical  force  density  pE 
is  nearly  equal  and  opposite  to  the  magnetic  one,  j  x  B/c  (as  expected 
from  the  force-free  condition) .  From  this  the  field  is  deduced  to  be 
nearly  curl-free  within  the  light  cylinder: 

rsinS  -  0(^-)  .  (4-50) 


As  mentioned,  near  the  light  cylinder  corotation  currents  become 
important  and  distort  an  assumed  dipole  field. 

Beyond  the  light  cylinder  the  force-free  condition  may  still  be 
assumed  -  the  restriction  that  8  <  1  will  then  allow  the  shape  of  the 
field  to  be  roughly  determined.  Equation  (4-43)  will  still  be  useful 
in  determining  the  charge  density,  but  in  (4-40)  the  value  of  K  cannot 
be  zero;  (4-40)  may  be  rewritten  as 
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Since  K  is  scalar  (cf.  eq. 


e<p  =  kb^ 


+ 


ursinQ 

c 


4-39)  Bp  =  KBp,  and 


n  ursinS 

2*.  -  — ~ 

BP  "  bp 


(4-51) 


(4-52) 


Assuming  0  <  8^  <  1,  some  manipulation  of  inequalities  yields 


so  that 


<  ursine  2 
Bp  ~  U  c  ; 


(4-53) 


.  .  .  .  ..  .  _  ursine  _ 

is  a  strict  inequality,  for  - - —  >  1. 

If  the  field  lines  remain  equipotentials ,  as  assumed,  they 
cannot  close  beyond  the  light  cylinder,  as  charges  could  then  corotate 
at  super luminal  speed.  Rather,  they  must  be  swept  back,  so  that 
| B^ /Bp |  grows  large.  Charged  particles  moving  on  such  open  lines  could 
escape  to  the  interstellar  medium. 

The  basic  structure  which  has  been  deduced  for  the  magnetic 
field  of  an  aligned  rotator  is  illustrated  in  Figure  9.  The  basic 
structure  of  the  corotating  magnetosphere,  under  the  assumptions  of 
negligible  space  charges  and  currents  and  a  force-free  condition,  have 
been  adopted  basically  unchanged  in  many  recent  models  (see  Manchester 
and  Taylor,  pp.  181-182).  However,  the  structure  of  the  polar  region 
and  of  the  wind  zone  (the  region  outside  the  CRM),  and  the  fitting 
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Figure  9.  Magnetic  field  structure  near 
a  pulsar  (aligned).  Poloidal  B  only  shown. 


together  of  these  structures  with  the  CRM,  vary  from  model  to  model. 

The  electrostatic  potential  for  a  rotating  magnetic  dipole  in 
vacuum  may  not  change  drastically  in  the  presence  of  a  thin  plasma. 

From  (4-5) ,  this  potential  at  the  surface  is 

2 

$(r  =  a , 6)  =  -  \  BQU)a  -  (3cos20  -  1),  (4-54) 

b  c 

and  thus  is  zero  for  cos0  =  3  2.  In  the  Goldreich-Julian  theory,  for 
>  0,  electrons  stream  out  from  the  surface  for  angles  smaller  than 
this,  and  protons  for  angles  larger.  The  field  line  based  at 
0  =  cos-1  (3-^)  is  called  the  critical  line.  This  is  reversed  for 
a) - i_i  <  0.  This  is,  however,  based  on  the  assumption  that  the  surface 
potential  for  the  vacuum  case  carries  over,  and  with  the  charge  density 
(4~44)  implies  that  in  certain  cases  charges  flow  through  regions  of 
opposite  space  charge  in  the  CRM.  There  is  also  the  problem  that 
(4-9)  and  (4-4)  imply  that  the  radial  force  on  a  charged  particle 
(constrained  to  move  on  a  magnetic  field  line)  would  be 

Fra  -  cos^Oq 

so  that  only  charges  of  q  <  0  (electrons)  would  be  accelerated  away 
from  the  surface.  The  fact  that  for  cos©  >  3-^  positive  charges  at  the 
pulsar  surface  have  larger  energy  than  they  would  if  allowed  to  escape 
to  infinity  does  not  matter,  the  local  field  will  not  allow  them  to 
stream  out  as  envisaged  by  Goldreich  and  Julian.  Again  if  u*p  <  0  the 
opposite  is  true  -  electrons  would  be  held  in  by  the  field.  The  role 
of  stellar  charge  in  removing  this  problem  will  be  discussed  in  later 


pages . 


> 


The  difference  in  ionic  and  electronic  work  functions, 
discussed  in  Chapter  3,  has  entered  some  recent  models  (Ruderman  and 
Sutherland  1975;  Jackson  1976)  in  an  important  role.  These  models  will 
be  discussed  in  the  following  pages.  The  ratio  of  ion  to  electron 
field  emission  currents  in  applied  electric  fields  of  the  same 
magnitude  will  very  roughly  be 


exp  (&* 


where  m^/me  and  <\>±/<be  are  the  mass  and  work  function  ratios,  and  both 
ion  and  electron  emission  are  assumed  to  follow  the  Fowler-Nordheim 
equation  (Gomer  1961,  p.  9).  The  ion  charge  is  Zc. 

Field  emission  of  electrons,  given  the  mass  and  work  function 
differences  from  ions,  should  be  very  easy  compared  to  that  of  ions. 
This  is  largely  due  to  the  tunneling  ability  of  the  less  massive 
electrons.  The  treatment  of  pulsar  surfaces  as  electron-emissive  but 
ion-retentive  is  justified  at  the  temperatures  which  characterize  most 
of  the  surface. 


Rarified  Magnetosphere 

Jackson  (1976)  considers  a  slowly  rotating  aligned  magnetic 
dipole  surrounded  by  a  tenuous  atmosphere  in  which  the  force-free 
condition  does  not  hold.  As  an  initial  approximation  the  vacuum 
electric  potential  is,  for  co*y  >  0, 


$(r  >  a,0) 


u)B0a3 

cr 


[e 


3  V 


2 

P2(cos0) ] , 


where  the  star  has  a  charge 


0)Boa 


3 


with  e  to  be  determined. 


(4-55) 


c 


The 
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magnetic  dipole  has  an  electric  monopole  superposed  on  it.  The  radial 
dependence  is  such  that  the  electric  monopole  dominates  the  motion  of 
charged  particles  at  large  radii. 

In  a  similar  manner  to  that  of  the  Goldreich-Julian  model, 
electrons  are  emitted  from  the  polar  region  if  0  <  e  <  1.  The  magnitude 
of  electron  emission  current  densities  will  be  much  greater  than  that 
for  ions  due  to  the  difference  in  work  functions.  At  the  surface 


a  2  o  a  5 

E*B  «  [e  -  (— )  cos20J (— )  cos0 


(4-56) 


so  that  positive  charges  are  emitted  for  0  >  cos-1(e^).  These  are 
bound  to  move  on  the  magnetic  surfaces  of  (4-45) ,  and  remain  within 


rmax  3  sin2 (cos'2 (  h))  1  -  e 


(4-57) 


of  the  star.  The  electrons,  emitted  near  the  pole,  can  move  to  very 
large  radii  where  they  feel  only  the  asymptotic  Coulomb  field  of  the 
star  until  they  merge  into  the  perfectly  conducting  interstellar  medium 
at  radius  R^.  In  the  steady  state  the  potentials  of  the  polar  cap  and 
that  of  the  interstellar  medium  have  been  equilibrated  by  the  electron 
current  so  that  the  charge  of  the  star  undergoes  no  further  changes. 
This  allows  e  to  be  determined  with  (4-55)  as 

_  %r  li 

e  -  —  [e  -  f] 
or 

e  “  i(1  -  4)_1  ~  1-  (4-58) 
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To  stress  the  zones  of  dominance  of  the  magnetic  dipole  and  electric 
monopole  terms,  respectively,  Jackson  defines  a  magnetic  flow  region 
(MFR) ,  where  |e  x  b|  >  E2  and  an  electric  flow  region  (EFR) ,  where 
|e  X  b|  <  E2  (see  also  Jackson  1975,  p.  583).  The  boundaries  of  these, 
and  the  surface  |e|  =  |b|,  are  approximately: 


MFR: 

EFR: 

E=B : 


2Ft  —  a  2 

sine  >  —  [1  -  e  H~ )  2P2  (cos0)  ] 
P£  r 

^3  (1  +  3cos20) 

r  >  —  p  0  - ; — - - 

2  *  sm9 

r  =  |  Pj,  (1  +  3cos20)^ 


(4-59) 


where  p£  =  -^  is  the  light  cylinder  radius. 

Although  once  e  has  taken  its  equilibrium  value  the  electrons 
will  no  longer  reach  the  interstellar  medium,  they  are  still  emitted 
at  the  pole.  Their  subsequent  motion  causes  current  loops  from  above 
the  pole  back  to  the  star.  Using  (4-55)  it  is  readily  shown  for  e  = 
that,  in  travelling  along  a  magnetic  field  line  near  the  star,  a 
particle  finds  itself  at  the  same  electrical  potential  as  at  the 
surface  (i.e.,  A$  =  0  when  it  has  reached  0  =  sin-1(y).  If,  however, 
it  enters  the  EFR  before  reaching  this  "standstill"  point,  it  will  feel 
an  electric  force  which  will  dominate  its  motion  and  cause  it  to  move 
to  magnetic  surfaces  with  a  larger  0O.  Jackson  then  envisages  motion 
along  these  surfaces  back  to  the  star  to  form  a  current  loop. 

As  a  higher-order  effect  of  rotation,  Jackson  considers  the 
atmospheric  charge  distribution.  As  (4-56)  shows  E*B  =  0  on  the  loci 

7T 

0  =  y  (the  equatorial  plane)  and 

r0  =  ae~^cos0. 


(4-60) 


u>|h* 


a  circle  of  radius  y-e  2  centred  at  Z  =  ye-'2.  Charges  near  one  of  these 
loci  would  feel  a  force  towards  these  E*B  =  0  surfaces.  If  any  energy 
dissipation  mechanism  is  available  the  particles  will  settle  onto  the 
E*B  =  0  surface  as  a  minimum  of  their  electrostatic  energy.  Jackson 
suggests  curvature  radiation  (to  be  discussed)  as  a  dissipation 
mechanism  to  allow  the  formation  of  an  electron  cloud  about  a  polar 
dome,  and  an  ion  cloud  in  the  equatorial  plane,  and  points  out  that  the 
exact  dynamics  of  the  particles  would  be  an  extremely  complicated 
problem.  This  picture  of  a  low  density  magnetosphere  is  given  in 
Figure  10.  The  negative  charges  along  E*B  =  0  form  a  dome,  above  and 
below  which  the  charge  density  should  be  very  small.  Since  Jackson 
considers  densities  which  are  undoubtedly  orders  of  magnitude  lower 
than  those  of  the  Goldreich-Julian  model,  vacuum  dipole  radiation  would 
be  the  main  energy  loss  mechanism  in  this  model  (see  Ostriker  and  Gunn 
1969;  Gunn  and  Ostriker  1970).  It  is  a  self-consistent  model  and 'is 
physically  reasonable  -  Jackson  suggests  that  the  basic  morphology  will 
remain  even  in  higher  density  regimes.  That  this  is  true  when  rotation 
is  rapid  and  the  light  cylinder  near  the  star  is  very  doubtful. 
Nevertheless  the  suggestion  that  the  stellar  charge  must  be  non-zero, 
which  also  characterizes  oblique  rotators  (Jackson  1978b),  and  that 
the  monopole  field  due  to  this  charge  should  dominate  motion  of  charged 
particles  at  large  radii,  is  likely  valid.  In  addition,  (4-56)  clearly 
defines  zones  of  electron  and  ion  emission  due  to  surface  fields, 
which  it  was  noted  is  a  feature  lacking  in  the  Goldreich-Julian  model. 


' 
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Figure  10.  Jackson  magnetosphere. 


TheRuderman- Sutherland  Model 


The  final  model  to  be  considered  here  in  some  detail  is  that  of 
Ruderman  and  Sutherland  (1975) ,  in  which  the  retention  of  ions  at  a 
pulsar  surface  shows  its  most  extreme  effects.  Pair  production  above 
polar  caps  is  a  major  feature  of  this  model,  as  it  is  of  the  model  of 
Sturrock  (1971). 

The  Goldreich-Julian  model's  general  features  as  described 
earlier  are  also  basic  to  this  model,  with  the  important  restriction 
that  wy  <  0,  so  that  the  poles  emit  ions  (to  the  extent  that  this 
occurs).  The  case  co*y  >  0,  the  "antipulsar,"  is  not  considered  here, 
since  the  role  of  ion  retention  is  then  not  so  clear. 

Most  of  the  open  lines  are  lines  onto  which  a  small  ion  current 
may  be  injected  and  thus  leave  the  star.  The  first  implication  is 
that  the  star  will  become  charged,  as  in  the  Jackson  model.  It  is  then 
not  clear  whether  positive  current  loops  from  the  pole  to  a  polar 
annulus,  or  a  current  stream  out  past  the  light  cylinder  would  form. 

In  either  case  a  region  where  due  either  to  a  strong  parallel  component 
of  the  electric  field  (Jackson) ,  or  the  need  to  resupply  to  the  magneto¬ 
sphere  charges  lost  through  the  light  cylinder  (Goldreich  and  Julian) , 
a  strong  positive  currently  radially  outward  from  the  polar  cap  exists, 
if  the  surface  supplies  ions.  Neglecting  the  effects  of  its  own 
magnetic  field,  the  current  flowing  from  the  polar  cap  can  have  a 
maximum  rate,  in  particles  per  second,  of 
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This  would  correspond  to  the  local  charge  density  of  (4-44)  moving 
outward  above  the  polar  cap  at  speed  c.  In  a  steady-state  situation 
where  the  surface  could  supply  ions  at  any  rate,  even  particle  fluxes 
as  high  as  that  in  (4-61)  could  flow  away  from  the  pole  without  the 
local  density  being  changed.  However,  if  ions  are  not  readily  avail¬ 
able,  the  region  above  the  polar  cap  will  be  depleted  of  particles 
completely.  This  is  due  to  the  electric  field  acting  on  the  charge- 
separated  plasma  above  the  polar  cap,  and  what  happens  to  particles 
cleared  out  of  the  gap  does  not  matter  at  this  stage  of  the  discussion. 

Ruderman  and  Sutherland  (1975)  have  solved  the  problem  of  a 
magnetosphere  separated  from  a  neutron  star  by  a  spherical  vacuum  gap 
of  height  h,  and  find  that  the  gap  has  a  field  and  total  potential  drop, 
at  the  pole,  of 


and 


(4-62) 


(4-63) 


respectively.  In  the  more  realistic  case  where  a  gap  occurs  only  at 
the  poles,  (4-62)  and  (4-63)  will  be  true  near  the  centre  of  the  gap, 
provided  h<<rp.  The  top  surface  of  the  gap  will  consist  of  a  plasma 
in  which  E*B  =  0.  In  addition,  this  condition  that  the  magnetic  field 
lines  be  equipotentials  is  met  on  the  side  of  the  gap,  so  that  as  the 
gap  grows  in  height  due  to  the  loss  of  positive  charges  from  the  region 
of  the  pole,  a  cylindrical  geometry,  with  this  condition  on  the  walls, 
becomes  appropriate  to  the  problem.  Using  this  geometry,  Ruderman  and 
Sutherland  find  that  a  maximum  potential  difference  of  -^fp 


can  be 


developed  in  a  gap. 

Michel  (1979)  discusses  gap  formation  as  a  basic  property  of 
pulsar  magnetospheres,  capable  of  removing  certain  basic  inconsistencies 
of  force-free  atmospheres.  The  Ruderman-Sutherland  model  is  a  special 
case  in  Michel's  discussion  although  in  his  paper  the  existence  of 
vacuum  gaps  is  postulated  without  establishing  a  physical  reason  for 
their  existence.  In  the  preceding  pages  the  ion  retention  of  the  pulsar 
surface  has  been  shown  to  allow  gap  formation,  so  that  Michel's  study 
of  the  structure  of  such  gaps  is  of  interest  in  that  it  provides  a 
complimentary  approach  to  the  problem. 

The  next  consideration  in  the  Ruderman-Sutherland  model  is  the 

breakdown  of  the  gap  by  pair  production  from  y-rays  of  energy  greater 

than  2mec2  moving  in  the  intense  magnetic  field.  It  is  pointed  out 

that  even  the  natural  background  flux  of  such  y-rays  exceeds  10 5  s  1 

onto  the  polar  cap.  As  pointed  out  by  Sturrock  (1971),  Erber  (1966) 

gave  equations  for  pair  production  in  intense  magnetic  fields.  The 

created  particles  move  in  the  gap  electric  field,  acquiring  an  energy, 

wBnh2e 

from  (4-63) ,  of  — ^ -  electron  volts,  which  for  the  canonical  pulsar 

parameters,  and  assuming  h  -  103  cm,  is  1011  eV  (Ruderman  and 
Sutherland  1975).  These  energetic  particles  in  turn  emit  curvature 
radiation  (see  below)  which  allows  the  cycle  to  start  again  and 
continue  to  build  up  a  large  spark  discharge. 

To  consider  in  further  detail  the  mechanism  of  the  spark,  we 
now  briefly  discuss  the  processes  involved.  In  the  intense  magnetic 
fields  which  are  found  near  the  neutron  star,  electrons  and  positrons 
will  predominantly  have  the  lowest  of  the  discreet  energies  given  by 


(3-12)  or  (3-14).  It  is  not  immediately  obvious  why  this  is  so, 
given  that  the  particle  energy  could  be  of  large  magnitude  even  before 
acceleration  by  the  field,  and  that  a  momentum  component  perpendicular 
to  the  field  could  be  present  and  of  large  magnitude  (Erber  1966) .  If 
the  energy  associated  with  this  perpendicular  (quantized)  motion  is  of 
the  same  order  as  that  with  motion  along  the  lines,  one  finds  that  if 
y2B  >>  6  x  1015  G,  radiation  formulae  for  strong  coupling  of  particle 
and  field  must  be  used  (Chiu  1972).  The  lifetime  of  an  electron  or 
positron  in  this  regime  has  an  upper  limit  of  order  10-19  (Bq/B)  s. 
Clearly  a  particle  will  leave  this  regime  rather  quickly,  entering  the 
quantum  regime  where  lifetimes  are  of  order  10“19  (Bq/B)2  s  (Chiu  1972), 
so  that  the  upper  levels  cannot  remain  populated.  The  observation  in 
Her  X-l  of  a  hard  X-ray  line  by  Trumper  and  his  colleagues  at  the  Max- 
Planck-Institut  (1977),  which  they  interpret  as  a  Landau-level 
transition  (see  Chapters  1  and  3)  could  be  due  to  this  decay  of  the 
perpendicular  energy  of  pair -produced  particles.  This  hypothesis  is 
advanced  cautiously  due  to  the  fact  that  gas  accretion  in  a  binary 
X-ray  pulsar  such  as  Her  X-l  (Manchester  and  Taylor  1977,  pp.  85-87) 
could  also  supply  particles  capable  of  producing  such  a  line.  The 
detection  of  such  lines  in  radio  pulsars  would  support  the  hypothesis 
that  pair  production  occurs,  since  it  would  provide  a  means  of 
populating  upper  levels.  What  is  clear,  however,  is  that  the  bulk  of 
the  particles  in  motion  along  field  lines  will  be  in  the  lowest  Landau 
level.  Their  radiation  can  no  longer  be  considered  classically,  and 
the  result  of  being  in  the  lowest  discreet  level  is  that  they  do  not 
radiate  due  to  their  motion  around  the  field  lines.  However,  in 


following  field  lines,  which  are  curved,  they  undergo  an  acceleration 
which  is  in  the  classical  regime  and  may  thus  be  treated  by  the  theory 
of  classical  synchrotron  radiation.  The  radius  of  curvature  entering 
the  synchrotron  formulae  will  be  pf,  the  local  radius  of  curvature  of 
the  field  line  which  the  particle  is  on.  The  characteristic  energy 
radiated  by  a  particle  in  such  motion  is  (Jackson  1975,  p.  675; 
Ruderman  and  Sutherland  1975) : 


Eph  =  -  2T3ftc/pf 


(4-64) 


and  they  are  emitted  within  an  angle  y_1  of  the  field  line.  Photons, 
either  from  external  sources  or  from  curvature  radiation  within  a 
region  of  intense  magnetic  field,  have  mean  free  path  due  to  pair 
production: 


for  Hu)  <  2mc2 


£  =  oo 


•ftw  >  2mc2  £  = 


A:.A._  JL  Bg.  exp  (A-) 

(e2/hc)  me  Bjl  p 


(4-65) 


(Ruderman  and  Sutherland  1975),  where  £  =  2 me"2  anc*  =  ®  sin6, 
where  0  is  the  angle  between  B  and  the  photon  propagation  vector. 

The  particle  multiplication  factor  needed  for  a  discharge  to 
occur  is  provided  by  the  many  photons  of  high  energy  emitted  by  a 
particle  as  curvature  radiation.  The  energy  loss  rate  (Jackson  1975, 
p.  661;  Ruderman  and  Sutherland  1975), 


P  =  Y4(t^) 


3  c3 


Pf 


(4-66) 


if  multiplied  by  he-1,  the  time  spent  traversing  the  gap,  gives  the 


energy  lost  by  one  particle  in  the  gap.  Dividing  by  the  mean  photon 
energy  from  (4-64)  the  number  of  photons  produced  is  approximated  by 

4  e2  h 

NPh  “  g  ^h  ^7  Y*  (4_67) 

(Note  error  in  eq.  25b,  Ruderman  and  Sutherland  1975.)  >>  1  is  the 

condition  for  a  discharge. 

The  gap  height  h  is  determined  by  requiring  that  1  =  h,  and  by 
using  the  mean  photon  energy  (4-64)  and  an  appropriate  Lorentz  factor 
corresponding  to  the  traversal  energy  of  the  gap,  producing  a 
transcendental  equation  for  h: 


where 


h2  =  1.45  exp  ah~7 


(4-68) 


1  =  1  /  eu)B  v  3  tic  1  1  _B_ 

a  8  me3  pf  2mc2  pf  Bq' 


(4-69) 


The  field  line  curvature  radius  pf  is  assumed  to  be  106  cm  in  the 
polar  gap,  typical  of  higher -order  multipole  fields  near  the  star.  The 
value  — B  is  that  of  Bj_  for  a  curvature  radiation  photon  after  going  a 
distance  h  from  its  point  of  emission  essentially  tangent  to  a  field 
line.  Ruderman  and  Sutherland  obtain  values  of  5”1  -  15  and 

3  Pf  2/7  3/7  -4/7 

h  =  5  x  103(qQg)  P  B12  cm  (4-70) 

for  the  gap.  Using  this  value  of  h,  the  photon  mean  energy  and  the  gap 
potential  may  be  expressed  in  more  fundamental  units: 


Eph  -  w3B3h6Pf  1  cc(pf5w3B  3)1/7 


(4-71) 


(4-72) 


1.6  x  1012B 
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Volts . 


For  a  canonical  pulsar  this  potential  difference  would  cause  electrons 
or  positrons  to  have  Lorentz  factors  of  order  106.  A  thin  gap  with 
this  potential  difference,  and  with  E*B  =  0  in  the  plasma  above  it  and 
at  the  top  edge,  will  have  a  parallel  electric  field  of  magnitude 

-  -  2o)B 

E-B  -  (h  -  z),  (4-73) 


where  z  is  the  height  in  the  gap  from  the  surface.  At  the  surface  this 
field  is  about  »  a  factor  h/a  -  10-3  smaller  than  (i.e.,  about 

10^  Vcm  !)  which  would  typify  a  vacuum  rotator.  Thus  with  the  binding 
energies  for  ions  discussed  in  Chapter  3,  ionic  field  emission  from 
polar  regions  beneath  such  vacuum  gaps  would  be  totally  negligible. 

For  the  parameters  just  obtained,  equation  (4-67)  gives  a 
multiplication  factor  of  50  photons/particle,  so  that  a  discharge  should 
easily  be  maintained. 

The  creation  of  pair -production  discharges  is  illustrated  in 
Figure  11,  where  a  major  problem  may  be  seen.  The  discharge  moves 
away  from  the  centre  of  curvature  of  the  field  lines,  and  in  the  case 
of  dipole  field  ends  up  on  axis.  Amplification  along  one  field  line 
does  not  occur,  so  that  a  shower,  instead  of  a  spark,  occurs.  Cheng 
and  Ruderman  (1977)  suggest  that  the  spark  picture  may  be  saved  by  the 
interaction  of  accelerated  electrons  with  the  surface,  which  will 
produce  a  y-ray  shower  leading  to  further  pe ir-creation  near  the  field 
line  on  which  the  original  pair  was  produced.  It  may  also  be  that  the 
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Figure  11.  Pair-production  discharges 
illustrating  drift  toward  magnetic  axis 
(after  Cheng  and  Ruderman  1977b) . 


synchrotron  radiation  which  may  accompany  the  birth  of  a  pair  (as 
discussed  earlier  in  this  chapter)  would  be  sufficiently  energetic 
while  the  particle  is  in  the  strong-coupling  regime  that  further  pairs 
could  be  produced  from  it.  It  is  worth  noting  that  if  y  -  10 5 
characterizes  the  energy  of  a  newly  formed  pair,  the  energy 
hwc  =  116  y2 ( Yq~3 )  keV  (Chiu  1971)  which  characterizes  the  critical 
frequency,  will  be  of  order  lOOmc2.  Erber  (1966)  notes  that  the  pair- 

g 

production  attenuation  coefficient  is  maximized  when  iiu)  =  12(-^-)mc2 
-  120mc2.  The  mean  free  path  for  such  photons  should  be  very  short 
and  thus  allow  further  pair  production  very  near  the  original  field 
line. 

The  gap  electric  field  accelerates  the  created  electrons 
toward  the  polar  cap,  and  the  positrons  toward  the  interstellar  medium. 
Fawley,  Arons,  and  Scharlemann  (1977)  claim  that  the  inward  electron 
flux  would  heat  the  surface  to  the  extent  that  free  ion  flow  (thermionic 
emission)  would  occur.  This  possibility,  and  the  results  of  the 
positron  flux  to  the  interstellar  medium,  will  now  be  considered. 

Ruderman  and  Sutherland  (1975)  give  as  an  estimate  of  the 
maximum  ion  flux  at  temperature  T: 

Pm  kT  ^  Eb 

F(T>  ■  v(v>  ex<>(-  (4-74) 

(see  also  Jenkins  and  Trodden  1965,  p.  71),  where  m  is  the  ion  mass 
and  pm  the  mass  density  at  the  surface.  Eg  is  the  ionic  binding  energy, 
essentially  the  cohesive  energy,  and  is  in  the  range  2-10  keV.  Due  to 
the  exponential  term,  only  when  kT  -  8.6  x  10  ®T  keV  also  is  of  this 
order  or  greater,  will  there  be  essentially  free  thermionic  ion 


emission.  The  failure  to  observe  X-ray  emission  has  ruled  out  pulsar 
surface  temperatures  of  this  order.  If  the  maximum  polar  current  of 
equation  (4-61)  were  due  to  thermionic  ion  emission,  a  temperature  of 
order  6  x  106  K  (Ruderman  and  Sutherland  1975)  would  be  required  at 
the  cap.  This  appears  to  be  observat ionally  ruled  out  by  an  absence  of 
X-ray  emission.  Of  course,  a  consistent  way  of  seeing  what  the 
temperatures  are,  is  to  calculate  them  based  on  the  model  being 
considered.  Ruderman  and  Sutherland  assume  inflow  of  electrons  at  the 
local  charge  density  of  a  force- free  atmosphere,  at  speed  c,  to  show 
what  temperatures  are  to  be  expected.  Assuming  blackbody  radiation  as 
the  only  loss  mechanism  for  the  ^AVpc  erg  cm-2s-1  flowing  into  the 
surface,  a  temperature  of  2.4  x  106  K  may  be  assigned  to  the  cap. 
Although  at  this  temperature  thermionic  ion  emission  could  not  supply 
ions  at  the  rate  of  equation  (4-61) ,  it  is  clear  that  under  these 
extreme  circumstances  ions  could  be  emitted,  thermionically ,  at  the 
polar  cap.  The  possibility  that  ion  emission  does  occur,  and  its 
influence  on  pair-production  discharges,  is  considered  in  Cheng  and 
Ruderman  (1977b).  However,  the  assumption  of  Fawley  et  al.  (1977)  that 
ion  emission  at  the  polar  caps  is  essentially  free  is  not  justified. 

This  is  seen  even  more  clearly  in  noting  that  a  pair -production 
discharge  is  needed  in  the  first  place  to  permit  the  heating  of  the  cap 
to  the  point  where  any  ions  are  emitted.  There  is  major  heating  only  in 
the  fastest  pulsars  (Cheng  and  Ruderman  1977c)  and  the  effects  of  the 
limited  ion  emission  in  others  will  not  be  discussed  here. 

An  observable  effect  of  positron  outflow  from  the  pulsar  (or  of 
Jackson's  return  currents)  should  be  the  presence  of  .511  MeV 


annihilation  radiation  from  pulsars.  Such  radiation  is  produced  by 
either  direct  annihilation  or  the  decay  of  singlet  positronium  which 
may  be  formed  by  the  positrons  with  electrons  from  the  interstellar 
medium.  Triplet  positronium,  which  decays  into  three  Y_rays>  giving  a 
continuum  spectrum,  could  be  formed  at  the  low  densities  in  the  inter¬ 
stellar  medium  (Ramaty  and  Lingenfelter  1979),  but  we  shall  assume  that 
it  is  absent. 

Again  let  us  assume,  as  an  extreme  case,  that  the  maximum  charge 
loss  of  equation  (4-61)  occurs  as  positrons  (there  would  also  be  a 
Goldreich-Julian  type  electron  wind  from  lower-latitude  field  lines  if 
a  Jackson  current  loop  is  not  present).  Every  positron  will  annihilate 
to  produce  two  photons  of  .511  MeV,  the  distribution  of  which  we  may 
assume  to  be  isotropic.  Then  at  earth,  a  distance  d  from  the  annihila¬ 
tion  region,  the  number  flux  of  .511  MeV  y-rays  will  be 

N  -  r  2  (4ud2)  1  cm~2s-1,  (4-75) 

P  e 

or  for  a  pulsar  with  rp  given  by  (4-47),  BQ  -  1012  G,  and  d  in  kpc: 

N  -  2  x  10“14p-3d-2  cm_2s-1.  (4—76) 

For  the  Crab  (P0531+21)  and  Vela  (P0833-45)  pulsars  this  would  imply 
that  about  10“ 10  photons  cm-2s_1  could  be  observed,  while  for  the 
nearby  pulsars  (d  =  .1  kpc)  P0950+08  and  P1929+10,  respectively 
1.2  x  10“ 10  and  1.7  x  10-10  photons  cm_2s"1  would  be  expected.  The 
present  limits  to  detection  at  .511  MeV  are  (liberally)  10“6  photons 
cm_2s_1  (HEAO-C  should  detect  10-4  to  10~5  cm“2s-1  over  the  range  .06 
to  10  MeV  (Lingenfelter  and  Ramaty  1978)).  The  upcoming  launches  of 


new  satellites  equipped  to  detect  y-ray  lines  should  allow  the  detection 
of  much  lower  levels  of  photon  flux  in  the  near  future.  The  four 
pulsars  mentioned  will  probably  be  the  easiest  to  detect  as  they 
maximize  P~3d“2. 

Time  Structure  of  Discharges 

As  the  discharge  proceeds,  the  gap  potential  is  reduced 
steadily.  The  time  scale  for  each  stage  of  pair-production  amplifica¬ 
tion  is  about  1/c  -  h/c  -  10*"3  s,  the  time  it  takes  a  curvature 
radiation  photon  to  move  through  its  mean  free  path.  This  is  the  time 
required  for  "tripling"  -  the  particle  which  emitted  the  curvature 
radiation,  plus  a  new  pair,  exist  after  this  time.  Once  enough  pairs 
have  been  produced  that  the  force-free  density  is  reached,  the  gap 
potential  will  be  zero.  Sometime  shortly  previous  to  this,  pair 
production,  and  thus  the  growth  of  the  discharge,  stop.  Ruderman  and 
Sutherland  (1975)  estimate  that  1015  pairs/cm2  must  be  produced,  which 
requires  about  35  (10 15  -  e35  -  335)  tripling  times,  or  10  ys.  This  is 
the  time  scale  of  the  micropulses  observed  in  many  pulsars  (see 
Chapter  1  and  Manchester  and  Taylor,  pp.  46-47).  This  time  scale  seems 
questionable:  if  each  particle  produces  50  photons,  not  a  tripling  but 

a  hundredfold  increase  occurs  at  each  step.  10  ys  is  an  upper  limit  to 
the  discharge  time,  if  this  simplistic  discharge  occurs  (ignoring 
currents) . 

Once  the  discharge  has  taken  place,  the  cycle  starts  over.  As 
this  restart  is  due  to  the  magnetosphere  demanding  particles  from  the 
polar  cap  either  for  acceleration  into  interstellar  space  or  for  a 


. 


Jackson-type  current,  it  is  important  to  note  that  this  demand  for 
charges  in  a  way  not  fully  explained  in  this  thesis  is  just  as 
fundamental  to  pair -production  through  polar  gap  formation  as  is  the 
basic  reluctance  of  the  surface  to  give  up  particles. 

It  is  clear  that  the  discharge  time  scale  would  be  much  less 
than  that  of  the  rotation  of  the  neutron  star,  so  that  many  discharges 
occur  in  one  period.  The  pulses  themselves  are  explained  by  a  beam 
above  the  magnetic  polar  cap(s)  swinging  by  us  as  the  (oblique)  rotator 
spins. 


A  further  timing  effect  of  this  theory  is  that  pulsars  should 
turn  off  once  their  period  becomes  too  long.  It  was  seen  (before 
equation  4-70)  that  £-1  <  15  is  a  condition  for  discharge.  From  this 
Ruderman  and  Sutherland  find  a  critical  period 


crit 
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(4-77) 


above  which  pair-production  discharges  do  not  occur.  Using  an  equation 
similar  to  (4-14),  and  with  a,  pf,  and  £  assumed  to  be  "typical,"  the 
result 

PP  <  const  (P/P)""13/5 


is  obtained.  In  1975  this  agreed  well  with  the  available  data,  and 
although  modifications  due  to  ion  emission  must  be  made  (Cheng  and 
Ruderman  1977b),  it  is  a  useful  and  testable  aspect  of  this  theory  that 
it  predicts  effects  which  should  be  seen  as  correlations  of  timing 


results. 


We  now  briefly  consider  the  oblique  rotator.  Although  Ruderman 
and  Sutherland  use  an  aligned  rotator  tilted  slightly  to  represent  an 
oblique  rotator  (and  in  so  doing  they  derive  expressions  for  the  drift 
periods  of  drifting  subpulses),  this  is  clearly  an  oversimplification. 

For  pair  production  to  occur  in  aligned  rotators  we  must  have  a  region 
of  ion  retention  above  which  positive  charges  are  quickly  swept  away  to 
form  a  gap.  In  the  case  where  coy  <  wy  <  0  and  with  no  stellar  charge, 
E*B  would  draw  electrons  from  the  star  between  the  rotation  and 
magnetic  equators,  and  ions  from  the  remainder  of  the  stellar  surface, 
but  if  the  star  has  charge  the  situation  becomes  much  more  complex 
(Jackson  1978).  Using  the  criterion  of  field  lines  closing  at  the 
light  cylinder  to  define  the  polar  cap  borders,  the  radius  of  the  polar 
cap  would  decrease  as  the  obliquity  increased,  then  increase  again  as 
X  tt/2.  Jackson  gives  the  value  of  the  parallel  component  of  E  at  the 
magnetic  pole  (again  wy  <  0)  as 

—  ~  B0wa 

E*B  - - cos0.  (4-79) 

c 

Since  for  0m  =  0,  cos0  =  cosx  and  is  constant  in  time,  the  situation  is 
similar  to  the  aligned  case,  but  the  polar  cap  may  not  be  spherical  due 
to  the  lack  of  symmetry,  and  the  field  becomes  small  for  large  x  (as 
an  extreme  case  consider  X  =  2”»  where  the  radial  electric  field 
precisely  at  the  magnetic  pole  is  zero,  and  the  ion-emissive  part  of 
the  polar  cap  is  a  semi-circle.  However  as  long  as  x  is  not  so  large 
that  the  parallel  electric  field  is  too  small,  pair  production  discharges 
would  still  be  expected  in  the  oblique  case. 


To  summarize,  an  aligned  magnetic  rotator  with  ion  retention 
would  be  expected  to  form  a  force-free  atmosphere  near  the  star,  and  to 
have  either  charged  particle  loss  along  open  field  lines  or  looping 
currents  in  the  outer  magnetosphere,  depending  on  the  extent  to  which 
it  is  charged.  In  the  special  case  that  rp  <  0  and  that  particle 
removal  from  the  polar  zone  is  very  efficient,  vacuum  gaps  will  form 
above  the  magnetic  poles,  and  discharge  through  pair  production.  The 
extension  to  oblique  rotators  is  plausible.  Observable  signs  that  this 
mechanism  is  present  would  be  lines  from  Landau  level  transitions  of 
newly  created  particles,  and  .511  MeV  radiation  at  a  level  below  current 
detection  limits.  The  radiation  at  other  wavelengths  which  such  a 
mechanism  produces  will  be  discussed  in  the  next  chapter. 


► 


CHAPTER  5 


PULSE  RADIATION  MECHANISMS 

This  chapter  will  mainly  extend  the  consideration  of  pair- 
production  discharges  to  show  that  they  could  lead  to  pulsar  radio 
emission. 

That  the  radio  emission  from  pulsars  is  coherent  may  be  deduced 

from  fluxes  observed,  combined  with  the  fastest  fluctuation  time 

scales,  of  order  100  ys  (Chapter  1).  The  size  associated  with  such 

fluctuations  is  of  order  10  km,  and  of  course  models  which  have  been 

discussed  are  also  of  about  this  size.  Knowing  the  distance  from 

dispersion  measures,  a  solid  angle  and  thus  intensity  Iv  may  be  assigned 

c2lv 

to  a  pulsar.  Radio  brightness  temperatures,  Tg  =  — which  are 
calculated,  are  typically  1023  to  1028  K  (Manchester  and  Taylor  1977, 
p.  201).  An  incoherent  radiation  mechanism  could  not  account  for 
brightness  temperatures  this  high.  Optical  and  X-ray  emission,  however, 
are  probably  incoherent,  with  brightness  temperatures  of  only  10 11  K. 

The  two  possibilities  for  production  of  coherent  radiation  are 
stimulated  emission  and  "antenna"  mechanisms.  Mechanisms  involving 
stimulated  emission  are  of  two  basic  types,  both  of  which  may  be  classed 
as  "maser  amplification":  direct  amplification  of  electromagnetic 
waves,  or  plasma  wave  amplification  followed  by  conversion  to  electro¬ 
magnetic  radiation  (Manchester  and  Taylor  1977,  p.  208).  Chiu  (1972; 
Chiu  and  Canuto  1971)  has  shown  that  the  first  mechanism  may  operate 
for  steady  electron  streaming  above  polar  caps,  due  to  a  population 


79 


inversion.  Soviet  workers  have  tended  to  favor  the  second  mechanism  and 
an  origin  of  the  pulsed  radiation  in  the  magnetosphere  relatively  far 
from  the  star  (a  "light  cylinder"  as  opposed  to  "polar  cap"  origin) 
(Ginzberg  and  Zheleznyakov  1975)  .  These  maser  mechanisms  generally  do 
not  require  special  spatial  distributions  of  radiating  particles  to  func¬ 
tion,  and  this  may  be  seen  as  an  asset.  However,  the  physical 
conditions  in  the  plasma  for  maser  processes  to  occur  are  somewhat 
restrictive  and  not  in  all  cases  consistent  with  the  generally  assumed 
magnetospheric  structure.  Light  cylinder  models  are  particularly  suspect 
due  to  possible  breakdown  of  corotation  (Mestel  1975;  Jackson  1976). 

The  aim  of  this  discussion  has  been  to  mention  that  from  maser  processes 
there  are  useful  emission  theories  which  should  be  considered  as  strong 
competitors  of  the  Ruderman-Sutherland  theory,  which  itself  is  an 
antenna  theory,  involving  the  formation  of  particle  bunches  of  diameters 
which  allow  their  curvature  radiation  to  add  coherently.  Fuller 
discussion  of  these  other  theories  is  to  be  found  in  Manchester  and 
Taylor's  book  (1977,  pp.  205-220). 

The  positron  beam  arising  from  pair  production  enters  the 
magnetosphere  above  the  gap,  still  moving  along  curved  magnetic  field 
lines  and  emitting  curvature  radiation.  The  production  of  pairs  as 
a  result  of  this  radiation  still  occurs  as  before,  but  since  in  the 
magnetosphere  E*B  =  0,  the  pairs  do  not  experience  acceleration. 

The  system  of  an  energetic  positron  beam  moving  through  a  plasma 
consisting  of  pair-produced  electrons  and  positrons  is  prone  to  the 
development  of  plasma  waves  by  the  two-stream  instability.  Cheng  and 
Ruderman  (1977a)  have  shown  that  the  electron-positron  plasma  must  stream 


in  order  to  maintain 


with  a  velocity  difference  of  (v+  -  v_)/c  -  ±  10"4 
the  charge  density  of  a  force-free  atmosphere.  They  find  that  the 
instability  which  develops  due  to  this  relative  streaming  then  has  a 
(laboratory)  frequency 


Re  a)  -  3  x  1013P“^B12  (^-)3/2  s"1 

and  growth  rate 

n  1  o  3/2 

Im  (i)  -  3  x  108  P^B12  (-)  s"1 


(5-1) 


(5-2) 


This  growth  rate  is  much  faster  than  that  associated  with  the  two-stream 

instability  due  to  the  energetic  positron  beam  moving  through  the  plasma 

(Ruderman  and  Sutherland  1975),  so  that  bunching  will  be  caused  by 

electron-positron  streaming.  This  is  due  to  the  fact  that  the  electron- 

"Y max  a  ^ 

positron  density  is  a  factor  —  (“)  >  50  more  than  that  of  the  primary 

beam,  due  to  the  large  number  of  curvature-radiation  photons  produced  by 
the  main  beam  of  Lorentz  factor  Ymax>  each  of  which  produces  pairs  with 
y+  =  •  The  fact  that  the  plasma  wave  instability  produces  bunching 

is  due  to  the  longitudinal  nature  of  E  in  such  a  wave  (Jackson  1975, 
p.  492).  This  bunching  by  separation  due  to  E  also  helps  prevent 
annihilations  from  occurring.  The  bunch  size  will  be  roughly  2ttco)“1, 
where  to  is  the  frequency  from  (5-1) .  This  size  is  of  the  order  of 
centimetres  so  that  coherent  radiation  will  be  produced  for  wavelengths 
similar  to  this.  Since  the  spectrum  of  curvature  radiation  rises  slowly 
as  (to/coc) 1  / 3  up  to  the  critical  frequency  toc,  but  falls  off  exponentially 
above  toc  (Jackson  1975,  pp.  676-677),  it  is  when  to  of  equation  (5-1)  is 
less  than  or  about  equal  to  toc  that  there  will  be  significant  radiation 


enhancement  by  coherence  (the  exponential  fall-off  of  the  spectrum 
overwhelms  the  coherence  increase,  which  is  only  a  power  law  with  index 
less  than  2  (see  e.g.  Jackson  1975,  p.  514)).  The  regions  where  this 
enhanced  emission  occurs  are  delineated  in  Figure  12.  Due  to  the 
dependence  on  r  of  equation  (5-1)  the  frequency  of  maximum  emission 
decreases  with  radius. 

At  a  given  frequency  (or  radius)  the  emission,  essentially 

tangent  to  the  field  lines,  occurs  between  angles  0m1'n  and  9_„v,  which 

iii-i-ii  max 

are,  in  degrees  (Manchester  and  Taylor  1977,  p.  223) 


min 


=  16P 


0.9 


w10 


1/3 


<W  =  16P-°-7  u)]0  1/3 


(5-3) 


where  wiq  is  the  radiation  frequency  in  units  of  1010  s_1 .  In  a  model 
where  the  beam  is  on  an  oblique  rotator  (subject  to  limitations  and 
distortions  discussed  in  Chapter  4)  the  sweep  through  the  line  of  sight 
will  determine  the  pulse  shape  (Figure  12b)  .  The  observed  variation 
of  pulse  width  with  period  seems  to  basically  agree  with  (5-3)  although 
the  frequency  variation  does  not  very  well  (Manchester  and  Taylor  1977, 
p.  223).  Spectra  and  polarization  of  the  pulses  are  also  predicted  by 
this  model  and  are  in  basic  but  not  very  precise  accord  with  observation 
(Ruderman  and  Sutherland  1975) . 

In  Chapter  4  the  expected  .511  MeV  annihilation  radiation  due  to 

positrons  (the  primary  beam)  escaping  to  the  interstellar  medium  was 

discussed.  We  have  shown  that  the  electron-positron  density  in  the 

^max  ,a 


magnetosphere  above  the  polar  caps  is  a  factor 


2tt_ 


(—t)  ~  50  greater 
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Figure  12.  (a)  Coherent  emission, 
regions  above  a  polar  cap 
(shaded) .  (from  Manchester  and 
Taylor  1977,  p.  222) 

(b)  Dependence  of  pulse  shape 
on  sweep  of  the  beam  across 


(b) 


line  of  sight 


than  that  of  the  primary  beam,  so  what  amount  of  yradiation  would  be 
expected  from  this  region?  Firstly,  it  should  be  noted  that  even  if 
the  annihilation  photons  completely  escape  the  magnetosphere,  the 
resulting  flux  of  say  fifty  times  the  estimates  of  Chapter  4  would  still 
be  below  current  detection  limits.  Secondly,  it  is  not  clear  that  the 
photons  will  escape  the  magnetosphere.  Hardee  (1977)  has  shown  that 
photons  of  energy  .511  MeV  (or  several  times  this)  should  have  a  mean 
free  path  length  allowing  them  easily  to  penetrate  the  entire  magneto¬ 
sphere,  considering  only  pair-production  attenuation.  In  the  highly 
relativistic  plasma,  collisions  may  be  energetic  enough  to  involve  pair 
annihilation  to  two  neutrinos  (instead  of  producing  yrays) ,  or  y-ray 
interactions  with  electrons  (or  positrons)  to  produce  neutrinos: 


e+  +  e“  -*■  v  +  v 

j.  +  — 

y  +  e~  -»•  e  +  v  +  v 


(5-5) 


(Chiu  1972)  .  Other  scattering  processes  may  also  be  present  to  prevent 
.511  MeV  line  radiation  from  the  magnetosphere,  but  have  not  been 
discussed  in  the  literature  in  this  context.  The  possible  pair 
annihilation  radiation  from  the  magnetosphere  leads  us  to  conclude  that 
the  photon  fluxes  calculated  in  Chapter  4  should  be  considered  minimum, 
rather  than  maximum,  estimates. 

The  pair-production  discharge  theory  has  been  shown  to  lead  to 
pulsar  radiation  characteristics  similar  to  those  observed,  in  addition 
to  being  consistent  with  the  structure  of  the  overall  inner  magneto¬ 
sphere  and  with  ion  retention.  The  ability  to  consistently  describe  a 
number  of  characteristics  of  pulsars  indicates  that  future  progress 


probably  depends  on  improving  the  understanding  of  certain  details  of 
this  theory. 
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APPENDIX  1 


Sources  of  Statistical  Data 

Manchester  and  Taylor  1977  (Appendix) 

Damashek,  Taylor,  and  Hulse  1978 
Manchester  et  al.  1978 

These  sources  present  data  for  321  pulsars.  Some  synthesis  of  these 
data  in  graphical  form  is  given  by  Manchester  (1979) . 
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APPENDIX  2 


Characteristic  Ages  of  20  Pulsars  and  Braking  Indices  of 
9  Pulsars  Based  on  Recent  Timing  Results 

Pulsar  frequencies  and  derivatives  were  observed  with  the 
Arecibo  radio  telescope  and  published  by  Gullahorn  and  Rankin  (1978b) . 

The  frequencies  and  first  derivatives  were  converted  to  period  and  period 
derivative  and  are  presented  as  Table  A-l,  along  with  the  characteristic 
age  in  millions  of  years  (Ma) .  These  data  are  also  presented  in 
Figures  5  and  6,  where  they  are  seen  to  blend  in  with  older  data 
without  introducing  any  significant  statistical  changes  to  those  data. 
Table  A-2  presents  frequencies  and  their  first  and  second  derivatives 
for  10  pulsars  -  the  first  nine  represent  new  results,  also  from  the 
(1978b)  paper  of  Gullahorn  and  Rankin.  The  last  listing,  given  for 
comparison,  is  the  Crab  Pulsar,  where  the  derivatives  are  well  known  due 
to  their  comparatively  large  magnitude  (Groth  1975) .  The  braking 
indices  listed  in  the  last  column  have  been  calculated  from  the  data  in 
the  table.  The  fact  that  they  seem  unreasonable  in  sign  and  magnitude 
(when  compared  to  the  Crab  or  to  theoretical  values  of  about  3)  casts 
considerable  doubt  on  the  second  derivatives  determined  by  Gullahorn 
and  Rankin.  These  authors  are  currently  preparing  a  series  of  papers 
discussing  these  results. 
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Table  A-l.  Periods,  Period  First  Derivatives,  and  Characteristic  Ages 
Tg  =  ^P/P  for  20  Pulsars,  Calculated  from  Data  by  Gullahorn  and  Rankin 
(1978b)  . 


p 

PCs) 

P(10"15ss"1) 

T6 (Ma) 

1900405 

.7465697 

12.896013 

.917 

1906409 

.8302699 

.0987758 

133.2 

1907412 

1.4417376 

8.2446963 

2.77 

1911413 

.5214722 

.8052111 

10.3 

1911411 

.6009975 

.6555952 

14.5 

19134167 

1.6162314 

.4076089 

62.9 

1913410 

.4045304 

15.244140 

.421 

1914409 

.2702529 

2.5182870 

1.70 

1914413 

.2818403 

3.6172045 

1.235 

1916414 

1.1808838 

211.35010 

.089 

1919+14 

.6181797 

5.6120633 

1.75 

1924+16 

.5798119 

18.003483 

.511 

1924+14 

1.3249220 

.2232364 

94.1 

1925+22 

1.4310666 

.7705831 

29.4 

1927+13 

.7600320 

3.6609526 

3.29 

1929+20 

.2682148 

4.1853985 

1.016 

1929+20 

.2682149 

4.1786406 

1.018 

1930+22 

.1444279 

57.779728 

.040 

1933+15 

.9673384 

4.0386196 

3.798 

1944+22 

1.3344500 

.8886812 

23.8 

2028+22 

.6305121 

.3837774 

11.3 

Table  A-2 


.  Frequencies ,  First  and  Second  Derivatives ,  and  Braking 
Indices  n  =  —■  for  9  Pulsars,  Calculated  from  Data  by  Gullahom  and 
Rankin  (1978b),  and  for  the  Crab  Pulsar  P0831+21,  from  Groth  (1975). 


p 

v  (Hz) 

v(10_15Hz  s_1) 

^(10"25Hz  s" 

2)  n 

0540+23 

4.0656 

-255.022 

3.953 

24.71 

0950408 

3.9516 

-3.57820 

-1.704 

-12664.70 

1541+09 

1.3361 

-.768157 

-1.116 

-252698.8 

1859403 

1.5257 

-17.42847 

1.100 

552.51 

1907+00 

.983337 

-5.332848 

-2.490 

-8609.61 

1915+13 

5.13805 

-190.1523 

-6.440 

-91.51 

1929+10 

4.41468 

-22.54446 

-3.18 

-2762.16 

2002+31 

.47366046 

-16.7302 

7.126 

1205.93 

2020+28 

2.912049 

-16.08157 

1.028 

1157.74 

0831+21 

30.2058647 

-3.855969  xlO5 

1.234084  xlO 

5  2.5071 

APPENDIX  3 


Very  Recent  Results 

Among  the  most  interesting  results  communicated  at  the  General 

Assembly  of  the  I.A.U.  in  Montreal  in  August  was  the  discovery,  among 

the  Molonglo  II  pulsars,  of  a  new  binary  pulsar.  In  contrast  to  the 

short-period  P1913+16,  the  new  system,  P0820+02  has  a  period  of 
ci  s  in  i  • 

1265  d,  and  - -  of  171  s  (if  P  =  0) ,  according  to  B.  J.  Robinson, 

so  that  the  easy  observation  of  relativistic  effects  is  not  expected. 
Nevertheless,  the  discovery  of  this  system  brings  the  "statistics"  of 
binary  radio  pulsars  (2  among  328  pulsars)  back  into  line  with  the 
expectations  of  binary  evolution  theory  (see  Chapter  2  and  reference 
therein).  J.  H.  Taylor  announced  that  the  pulse  shape  changes 
"observed"  in  P1913+16  (Taylor  et  al.  1978)  are  probably  an  observa¬ 
tional  effect.  He  also  advanced  the  idea  that  interpulses  may  not 
represent  emission  from  two  polar  caps  but  simply  be  an  extreme  case  of 
wide  pulses  -  a  pulse-width  histogram  supports  this  hypothesis. 

A.  G.  Lyne  presented  Jodrell  Bank  astrometry  results  on  9 
pulsars.  P1920+10  is  located  69  ±  26  pc  away,  but  for  the  8  others 
studied,  the  errors  are  larger  than  the  measurements  (the  parallax  is 
of  order  10“3  arcsec  and  is  measured  through  interferometry  with  a 
127  km  baseline) .  Distances  based  on  dispersion  measure  are  comparable 
to  those  of  these  parallaxes  if  a  value  of  <Ne>  =  .03  cm-3  is  adopted 
for  the  electron  density  (see  also  the  letter  by  Salter,  Lyne,  and 
Anderson  in  Nature,  vol.  280,  pp.  477-478). 


104 


The  theoretical  paper  of  most  relevance  to  the  Ruderman- 

Sutherland  class  of  models  was  presented  by  V.  Radhakrishnan.  Assuming 

that  pair-production  discharges  in  a  polar  cap  gap  are  initiated  solely 

by  background  y-rays,  and  noting  that  the  path  length  (see  equation 

4-65)  is  a  decreasing  function  of  the  photon  energy,  he  finds  that  the 

magnetosphere  acts  as  a  low-pass  filter,  on  photons.  The  most 

energetic  of  these  photons  entering  the  gap  will  have  the  most  chance 

of  initiating  discharges,  so  that  in  a  given  magnetic  field  a  band  of 

photon  energies  is  selected  as  possible  initiators  of  discharges. 

Conversely,  given  the  natural  background  spectrum,  one  can  find  that 

fields  of  Bj2  =  1  to  3  are  needed  for  pair-production.  This  correlates 

•  r. 

well  with  fields  determined  from  (PP) 2. 

A  final  recent  report  of  major  importance  is  the  "Possible 
optical  observation  of  the  companion  star  to  the  binary  pulsar 
PSR1913+16"  (Nature  280,  pp.  367-370),  in  which  P.  Crane,  J.  Nelson, 
and  J.  Tyson  report  observing  a  star  of  Mr  =  20.9  at  the  precise 
coordinates  of  the  pulsar.  Their  interpretation  of  this  object  as  a 
helium  (i.e.,  helium-burning  main-sequence)  star  would  introduce  tidal 
effects  to  the  binary  system  which  disallow  the  interpretation  (Taylor 
et  al.  1978)  of  timing  effects  as  due  to  gravitational  radiation.  The 
suggested  confirmation  (spectroscopically)  of  this  object  as  a  helium 
star  companion  to  PSR1913+16  is  of  obvious  importance. 


